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ABSTRACT 


The theory of varied flow was first postulated in a complete 
id comprehensive manner by J. M. Belanger. The fundamental 
ainciples were so well covered that little has since been added 
) modify the differential equation of gradually varied flow in 
S original form. All subsequent investigations have been con- 
ed chiefly to the methods of solution of the differential 
uation by direct, approximate, and graphical integration. 

In this bulletin three well-known methods of integration 
e discussed with special reference to their advantages and 
mitations, and a new method is introduced. Recent investiga- 
ons have tended to favor direct integration in order to dis- 
nse with the cumbersome succession of steps required by the 
proximate and graphical methods. Following the current trend, 

ore convenient procedure is introduced in which limitations on 
€ existing methods of direct integration are removed, but 
bodying at the same time all their advantages. The practical 
plication of the proposed method to the determination of 
rface profiles of flow in various channels is given. In the 
ecial cases of the:rectangular and circular channels charts 

e been prepared to facilitate the profile computations. 

The validity of the proposed method is determined by comparing 
e theoretical results with actual profile observation under 

ious conditions of flow in experimental flumes. The compari- 

of test data is also extended to the theoretical results 
puted by the more widely accepted methods. 
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1. INTRODUCTION 


The study of the hydraulics of flow in open channels is not 
jlways subject to analysis in an exact or rigorous manner, 
ecause of the large number of fundamental variables involved. 


jermits the use of models in the analysis of flow characteristics 
f prototype hydraulic structures. That the results obtained 
rom such simplification cannot be regarded as better than ap- 
roximate is evident, but frequently the limitations imposed by 
his simplification are compatible with the desired accuracy of 
fesults in practice. 

In the case of steady flow in open channels the forces present 
re those due to gravity, inertia, viscosity, elasticity, and 
urface tension, the last two being inconsequential except in 
mall models. Gravity forces are generally taken into considera- 
ion, while on the other hand either the viscosity or inertia 
orces may be neglected in the solution of many practical prob- 
ems in open-channel flow. These problems involve either one of 
wo broad classes of motion: rapidly varied flow in which gravity 
nd inertia are the predominant forces; and gradually varied 
ow, in which the flow pattern is determined primarily by grav- 
ty and viscosity forces. Rapidly varied flow is usually associ- 
ted with rapid or abrupt changes in the depth of flow occurring 
nh a comparatively short distance, such as the hydraulic drop or 
he hydraulic jump. Gradually varied flow, on the other hand, is 
haracterized by the absence of a disrupted surface over a con- 
iderable reach of channel, the changes of depth being so gradual 
hat the flow may be considered essentially parallel. Most of 
he problems in gradually varied flow are concerned with the 
etermination of surface profiles of flow in both natural and 
niform channels. The present work is confined mainly to this 
hase of open-channel flow. 


. Purpose and Scope of Research 


The purpose of this research was to develop a method of 
omputing surface profiles of gradually varied flow in uniform 
hannels whereby the results could be obtained more directly 
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than are possible with the existing methods. To determine tl 
accuracy of the proposed method, profiles of flow were recorde 
in experimental channels of rectangular, trapezoidal, and circ 
lar cross-sections under various conditions of flow. The tes 
data are compared with the results computed by the propose 
method ands also by other methods that are in use. This investi 
gation was confined mainly to backwater and downdrop curves 1 
uniform channels on mild slopes, inasmuch as they form the mor 
important classes of surface profiles encountered in practice. 


2. Notation 


The nomenclature used in this bulletin conforms essential] 
to that prepared by the Committee of the American Standards Asse 
ciation on Letter Symbols for Hydraulics, approved in January 
1942. The symbols, defined when they first appear, are summé 
rized below for reference. 


A: total cross-section area 
b: width of channel 
b6.: width of water surface 


C,Cy: coefficient of flow (Chezy) 
Cy: exponent in Mononobe’s functions 


C_: constant in the proposed varied-flow 
equation to account for nonuniform 
velocity distribution 


C,: a parameter in the proposed varied- 
flow equation 


D: diameter of circular conduit 
E: energy per unit weight 
E : mean energy per unit weight 
E.: specific energy 
G: os (A_)2 b t 
iy Renan 
gravitational acceleration 
potential head = y + Z 


conveyance of a cross-section 


es PY os 0 


a constant in the proposed varied-flow 
equation to account for eddy losses 


l: exponent in the wetted perimeter-depth 
relationship 
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M: exponent in the area-depth relationship 
N: hydraulic exponent 


n: coefficient of roughness (Kutter and 
Manning) 


»: Wetted perimeter of cross-sectional area 
Q: discharge rate of flow 
R: hydraulic radius of cross-sectional area 


r: ratio of depths, y/o 

age BLS ye 17 5 

S: slope of energy grade line 
Sat. Critieal (slope 


friction slope 


i 
Sy: slope of channel bed, tan ay 
V 


mean velocity, Q/A 
x: horizontal length 
y: vertical depth from bottom of channel 
y-: Belanger critical depth 
Yo: depth for uniform flow 
Z: elevation above datum 


ay: angle of an inclined channel bed 
measured from the horizontal 


B: Bakhmeteff varied-flow function 


A: used with x, y, r, B to indicate the 
difference between two particular 
values 


g: function of 
@,: Bresse backwater function 
B 
ee ae Mononobe backwater function 
Ya» Wa: Mononobe dropdown function 
Except for the classification of surface profiles, the sub- 
cripts ] and 2 are used to refer, respectively, to upstream and 


ownstream sections; critical conditions are designated by the 
ubscript c, and uniform conditions by the subscript 0. 


11. STEADY FLOW IN OPEN CHANNELS 


3. Uniform Flow and Varied Flow 


Flow in a channel is called uniform when certain elements 
related to the flow and the channel remain constant from sectior: 
to section. These elements are the depth of flow, the cross- 
sectional area of flow, the velocity of {low and the slope of the 
channel. Strictly speaking then, uniform flow can take place 
only in prismatic channels of constant slope. Although variec 
or nonuniform flow is the rule in nature, it can also be made tc 
take place in prismatic channels. 

This bulletin deals primarily with certain phenomena connectec 
with varied flow. Varied flow must also be distinguished fron 
variable flow. In varied flow the velocity and depth, while 
varying from section to section, do not vary with time. In 
variable flow, however, these elements vary with time as in the 
case of waves. There are two main types of varied flow - one if 
which the depth of the water increases gradually downstream, anc 
the other in which the depth of the water decreases gradually 
downstream. The frst type gives rise to a surface profile curve 
which is called a backwater curve. The second type results in € 
dropdown curve. 

It will sometimes happen that the conditions of flow are suck 
that there may be a sudden break in the surface of the backwater 
curve. We then get what is called a "hydraulic jump." Similarly 
a sudden break in the surface of a dropdown curve gives rise 
to a "hydraulic drop." 


4. Energy Relations in Varied Flow 


_ In order to derive the equation of gradually varied flow, it 
is first necessary to review the various energy relations exist- 
ing in uniform channels during the gradually varied flow. Con- 
sider a liquid flowing in a prismatic channel at an angle a») and 
having a constant slope S, = tan ay (Fig. 1). Let it be assumed 
that the velocity at any particular section is uniform over the 
cross section. Then, at any cross section, if E, equals the mean 
energy per unit weight of liquid, | 


er aay hae ees i 
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Where h = y + Z and V2/2g represents the average kinetic energy 
of the stream. 


Differentiating with respect to x, the displacement coordinate, 
Ge ecah da Ve 
Ae GO) (2) 
dx dx die 28 


Consider the flow to be uniform; i.e., V is constant. Then Eq. 2 
becomes 


dE dh 


Lie cpp 
ax dx ot 


Equation 3 states that in uniform flow, the rate of change of 
energy is equal to the rate of change of the surface above the 


4 


> 
Nie =e 


: Fig. ia Surface Profile of Gradually Varied Flow 


datum. It should be remembered that dh/dx will be numerically 
negative, since the change of energy in the downstream direction 


dh 
ms negative. , Since in uniform flow, GZ = -So, the slope of the 
channel, 
dE 
Ss “So (4) 
dx 


From Eq. 4 it may be seen that in uniform flow the work done by 
gravity is just sufficient to overcome the resistances of the 
channel. 

In uniform flow the liquid flows at a constant depth for a 
given rate of discharge Q. This "normal depth" is designated by 
the symbol y,. It should be remembered that for a given channel 
and a given discharge, y, is one of the characteristics of the 
channel. 

Returning to Eq. 2 for the general case of nonuniform flow, V 
will be a function of x. At any depth of flow 
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where A is a function of y and Q is constant. Differentiating 
with respect to x 


d 
— (v4) = (5) 
ax 
Fig. 2. Evaluation of dA/dy 
Substituting Fq. 5 into Eq. 2 gives 
dE, dh Q? dA dy 
dx dx gA® dy dx (6) 


To evaluate dA/dy, consider Fig. 2, the cross-section of ¢é 
prismatic channel with a survace width 6,. Clearly, 


dA = b. dy or —e=b6 (7) 


dx dx MeA> dx (8) 


Next, dh/dx can be evaluated in terms of y and x if in Fig. @ 
the slope is considered so flat that sin a, can be assumed equal 
t 

to the slope S,) and cos a, can be assumed unity, then 


dh = Spdx - dy 


‘dx 4 DOS vars (9) 


Noting that dh/dx is negative, if Eq. 9 is substituted into Eq. & 
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dx 0 La we gA3 dx (10) 


yt is now necessary to find an expression for the resistance loss 
mn varied flow dE/dx, in order to solve Eq. 10 


. Evaluation of Resistance Losses in Varied Flow 


While it is probably true that the resistance losses of varied 
ow follow a somewhat different law than that of uniform flow, 
yn most practical cases where the change in depth is gradual, it 
s close to the truth if the following assumption is made: the 
ate of resistance losses in varied flow at any section of depth 
ris identical with the rate of loss which would take place if 


Fig. 3. Evaluation of dh/dx 


the flow in the same channel were uniform with the same discharge 
and at the same depth. Uniform flow under these conditions can 
take place only in a channel whose bottom slope is different from 
the bottom slope of the channel in which the gradually varied 
ow is occurring. The bottom slope of the uniform flow channel 
ay will give the desired rate of loss due to frictional resis- 
tances can be obtained from any empirical equation of uniform 
flow. This resistance slope S, can be obtained, for example, 


from the Manning formula'*’ in which the foot is used as the unit 
length. 
eee ea Oe 
Sp = gee R278’ * “tage aRt73 (11a) 


In this case the frictional losses S, are inversely proportional 
to some function of the depth of flow. For developmental pur- 
poses Eq. lla is used in the form 
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4 
ue 
f K2 


where K is called the "conveyance" of the channel and is a fune 
tion of depth. The meaning of Eq. 11b should be clearly under 
stood. It states that for a given channel with a given discharg 
there is a certain slope S, which is just sufficient to overcon 
the resistances. The bottom slope Sy) of the gradually varie 
flow channel may be greater or less than S,,. The frictional re 
sistance S; is numerically equal to dE/dx but is different 1 
Sign, since the former is a rate of loss of energy whereas tlk 
latter is expressed as a rate of change of energy. Using tt 
relationships expressed in Eqs. 10 and 11, the equation of grac 
ually varied flow in uniform channels can be formulated. 


(11b) 


11. EQUATION OF GRADUALLY VARIED FLOW 
IN UNIFORM CHANNELS 


ss Formulation of the Equation of Varied Flow 
Substituting “S from Eq. 11 for dE/dx in Eq. 10 gives 


d y d 
Q iy 


SS ee Shute, 
i 0 ¥ dx w gA? dx (12a) 


r transposing 


a Routes 


2 
dx 6 Q 


(12b) 


dx = ——————- dy (13) 


iquation 13 expresses the change in x in terms of the change in 
* and functions of y and is the required differential equation 
or nonuniform flow in prismatic channels. Although Eq. 13 is a 
rst-order differential equation with variables x and y sepa- 
ated, its solution cannot be obtained by quadrature. For this 
€ason various simplifying assumptions have to be made, and dif- 
erent authorities have suggested different assumptions, more or 
ess useful from a practical standpoint. Before describing these 
olutions it is necessary to consider briefly the limitations 


rf Eq.) 13. 
| 


. Limitations of Applicability of the Varied-Flow Equation 
The equation of gradually varied flow was derived from certain 

implifying assumptions; therefore, it may not be valid beyond 
he premises upon which the derivation is based. These assump- 
ons are as follows: 

(a) Steady flow in uniform channels on flat slopes 

(b) Hydrostatic distribution of pressure intensity 

over a normal section 
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(c) Uniform distribution of velocity across the 
entire section of flow. 


The applicability of Eq. 13 is necessarily limited to fle 
slopes because in the derivation'*’, sin a) is assumed to | 
equal to tan ao (Sy) and cos oy and cos2 a) are assumed to | 


unity. In other words, the depth of flow was assumed to be tk 
same whether the vertical or normal direction was used. Similan 
ly, the distance along the bottom of the inclined channel we 
assumed to be the same as its horizontal projection. For 

channel slope of 1:10 the value of cos? ag is 0.99 resulting i 
an error of only 1 percent. 

Strictly speaking, hydrostatic distribution of pressuy 
intensity is attained only when the normal component of accel 
eration is zero. This condition imposes the restrictions thé 
the stream lines be neither divergent nor curvilinear'?). Fe 
gradually varied flow in uniform channels, these conditions a1 
essentially met since the changes in depth are so gradual the 
the flow may be considered substantially parallel. If the fle 
involves rapid changes in depth or direction, the acceleratic 
effects cannot be neglected and the equation of gradually varie 
flow is not applicable. 

The assumption of uniform velocity distribution general 
involves little error, particularly in cases where the velocii 
head is small compared with the depth, as is the case in graa 
ually varied flow on mild slopes. This is due to the fact thé 
the velocity-correction factor is actually very close to unity™ 
Values of this factor under various conditions have been found 
vary between 1.01 and 1.12, averaging about 1.05. If the velow 
ity distribution deviates substantially from unity, then tk 
velocity-correction factor must be considered in the derivatic 
of the gradually varied-flow equation. 


8. Specific Energy, Critical Flow, and Critical Slope 
In order to obtain a better concept of the relationshi 
between S,) and Sf Bernoulli’s equation for expressing tk 


conditions between two cross-sections of a stream, 1] and { 
may be stated in the form | 


8 Sel ee We ot (K.E.), = UETES ) » + (K.E.), +(losses))_, (14) 


In the equation, P.E. stands for the potential energy head, | 
K.E.. for the kinetic energy head, V2/2g; and the lossés ated 
sumed to be equal to the ficient ab head losses between section 
l and 2. Transposing 
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(PSE ~) 5 - (P.E.)5 - (losses) ,_, = (K.E.), - CRE) (i) 


or a length of channel dx, we may write this relationship as 


Soy ase 3 S; GeRe INK Beh g um (KE) 5 (16) 
quation 16 indicates that if S) is greater than S,, the work 
one by gravity is greater than that required to overcome fric- 
ion and that some energy will be added to the stream. If, how- 
ver, Sy 1s greater than S), some energy will be withdrawn from 
he store of energy in the stream. The former condition will be 
epresented by a rise of the water surface downstream and the 
atter by a fall of the water surface downstream. 

These relations can be made clearer by the concept of "speci- 
c energy" of the flow, E,. If we refer to the bed of the stream 
s datum, the specific energy of the stream at any section is: 


Vas 
Ey i + 2g (17) 


here y is the depth of flow at the section. The distinction 
etween Eqs. 1 and 17 is important. Equation ] is a measure of 
he mean energy of the stream above any datum; Eq. 17 is a 
easure of the specific energy of the stream where the bottom 
f the channel is used as the datum. In uniform flow, E, 1s con- 
tant. In varied flow, E, will change from section to section. 
For dE/dx positive, the specific energy increases downstream. 
1s means that the work of gravity exceeds that required to 
vercome resistance losses, and the excess adds to the specific 
nergy. If the work of gravity is insufficient to overcome the 
frictional resistances, part of the specific energy 1s used to 
ake up the difference. 

_ By plotting E, (Eq. 17) against y it may be shown that for any 
channel with a given discharge, there is a depth y for which E, 
is a minimum. This depth is called the critical depth, y,. A 
given discharge Q may flow in a channel at various depths, y. 
For each value of y there corresponds a value of the specific 
energy given by Eq. 17. But under no circumstances will the 
content of specific energy fall below a minimum value of E,, and 
this minimum is attained at the critical depth, y,. Substituting 
JA for V in Eq. 17, differentiating with respect to y, and 
equating to zero gives 
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w 


Since 24 = 6 (Baa 
dy 


dE b Q2 
peek — i pre = 0 (18) 
dy gA 
or 3 9 
A 
ME (19) 
6 g 


In other words, the depth of flow for a given discharge is critil 
cal for the conditions expressed in Eq. 19. For a given diss 
charge, the slope that will make the flow critical at the norma 
depth is called the critical slope and is designated S_. 


9. Classification of Surface Profiles 


Longitudinal profiles of flow in open channels may be conven: 
iently classified according to the bed slope Sj and the depth og 
flow y in the following manner. !?’ 


Table | 
Types of Surface Profiles in Uniform Channels 


Slope Designation Relation of Depth, y, 
to Uniform and 
Critical Depths 


Adverse, Ay Ose Dyas a 
So < 0 Ag JoSiyony 
Horizontal, Ho Yours Yeas 
Mild, My VSTIG HH 
Sale olde ne Y> > Ye q 
M3 SQ ata 0 
Critical, Cy 1 > Yo = Y0 
So = S,> 0 C3 YO.= %o> ¥ 
Steep S} YS eI 
Sp > S, > 0 So Mer. Tig 


S3 Ae te eae | 
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If So < 0, the channel is termed adverse A; if Oj tes kt 
s horizontal H; and if S) > 0, it is called sustaining and is 
urther classified by the degree of slope. Uniform motion in 
pen channels of sustaining slopes may proceed in the tranquil, 
Pitical, or rapid state depending on whether the uniform depth 
f flow Yo 1S greater than, equal to, or less than the critical 
epth, y,. Sustaining slopes producing uniform motion in the 
ranquil, critical, and rapid states are said to be, respectively 
ild M, critical C, and steep S. 


Fig. 4 Longitudinal Profiles of M, and Mo Curves 


If the actual depth of flow y lies above both the uniform- 
depth and critical-depth lines drawn parallel to the bottom of 
the channel, it is of type 1; if between these lines it is of 
type 2; and if below both lines, it is of type 3. The twelve 
types of surface profile are outlined in Table I. 

Of these surface profiles, the more common types are those on 
ild slopes M, and M, (See Fig. 4). The former refers to the 
portion of the water surface lying above the uniform depth, and 
the latter to that lying between the uniform and critical depths. 
The M, curve, commonly known as the "backwater" curve, is pro- 
duced when the downstream depth of flow is raised above the uni- 
form depth by a dam or other obstruction. The M,, or "dropdown" 
curve, results when the downstream depth is lowered below the 
uniform depth. The upstream and downstream limits of the M, and 
Mo curves are reached when y = yy ees andy = Yor. a 
respectively. The conditions of flow existing at these limits 
nay be found by consideration of the gradually varied-flow equa- 
mon. That is: 
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(12b) 


eA? 


At the upstream limit, y = yg; the friction slope S, and than 
of the channel Sy are identical; therefore, the left side o: 
Eq. 12b becomes zero. This condition indicates that both the MJ 
and M, curves approach this limit asymptotically. Theoretically/ 
the length of these curves extends to infinity. For practicai 
purposes, however, these curves may be assumed to end at a depth 
slightly above (for the My curve) or below (for the M, curve) tha 
uniform depth. 

At the downstream limit for the M, curve, y = %, the friction 
slope S, becomes zero, while the denominator of the right side oi 
Eq. 12b has a magnitude of unity. The left side of the equation, 
consequently, is equal to S,. Thus, the M, curve approaches tha 
downstream limit horizontally. 

At the downstream limit of the M, curve, y = y,, the denonina~ 
tor of the right side of the equation becomes zero, so that the 
left side approaches negative infinity as a limit. Consequently, 
the M, curve should end downward in a vertical direction. 

The shapes of the M,- and M,-types of surface profile are 
shown in Fig. 4 in greatly reduced horizontal scale. It shoule 
be emphasized that the limiting conditions serve only the useful 
purpose of establishing the general shape of the different types 
of surface curves. Before many of the limits are reached, the 
assumption of parallel flow may no longer be valid. Equation 13, 
then, cannot be correctly applied to investigate the existing 
state of flow under those conditions. 


V. INTEGRATION OF THE GRADUALLY VARIED-FLOW EQUATION: 
EXISTING METHODS 


: General Review 


The theory of varied flow was first postulated in a complete 
nd comprehensive manner by J. N. Belanger “°) His paper, pub- 
ished in 1828, "Essai sur la solution numerique de quelques 
roblemes, relatives au movement permanent des eaux courantes, " 
Jontains, in part, the general differential equation for gradual- 
y varied flow and a method of solution by successive approxima- 
1ons. The fundamental principles were so well covered that 
ittle has since been added to modify the differential equation 
mh its original form. All subsequent investigations have been 
jonfined mainly to the methods of solution of the gradually varied 
ow equation by direct, approximate, and graphical integration. 

The earlier attempts to obtain an analytical solution by 
irect integration of the differential equation were restricted 
o flow in channels of special form. The case of a rectangular 
hannel of great width was treated by Dupuit (1848) and in a 
omewhat different manner by Ruhlmann (1880), both ignoring the 
flect due to changes of velocity.'®’ The same case was pre- 
ented in complete form by Bresse (1860) and subsequently by 
rashof.'’’ These methods of integration were based on the Chezy 
quation with a constant coefficient in evaluating the friction 
lope. More recent solutions, based on a variable coefficient in 
he Chezy equation, were due to Schaffernak (1913), Ehrenberger 
1914), Baticle (1921), Kozeny (1928), Schoklitsch (1920), ‘’’ and 
under (1943). '° The case of a broad parabolic channel was 
reated by Tolkmitt (1898).‘’’ The first attempt to arrive at a 
olution suitable for any type of cross-section was begun by 
akhmetefl (1932), ‘9’ and subsequently presented in a more 
omplete form by Mononobe (1936). ‘7’ 

In the solution of the differential equation for gradually 
aried flow by approximate integration, the channel is first 
ivided into several reaches and the depth at the end of each 
“each determined by a trial-and-error procedure. A more direct 
solution was given by Husted (1924), °°’ in which the lengths of 
successive reaches were found for assumed increments of depth. 
Since the results are obtained by a series of successive computa- 
ions, this method is generally known as the "step method." 
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The gradually varied-flow equation may also be solved by 
graphical integration. ‘17"*?) The graphical methods involve 
plotting suitable variables in rectangular coordinates in such ; 
manner that the distance along the channel is given by the arez 
under the curve. The area is found by means of a planimeter. 

All three methods, direct, approximate and graphical integra: 
tion, are used to a varying extent in computing surface profiles 
of flow in uniform channels. 


ll. Integration of the Equation of Varied Flow 


Several methods, more or less exact, are available for tha 
integration of Eq. 13. Only the more important ones are de- 
scribed briefly here. These are the methods of graphical integrai 
tion, the step method, and the direct integration methods 0: 
Bakhmeteff and Mononobe. For discussing and comparing theses 
methods it 1S convenient to express Eq. 13 in a different form, 
From Eq. 13 


502 
ety set eg 
gA 
pital ee SE Oy 
So - Sy 
b Q? 
1 - w 
aie ap li thts <5 (20) 
1 pas a 
So 
But 
Q2 
Sy * AGIRI73 
2 
Peay ye . 4/3 
0 0 


where R is the hydraulic radius at normal depth y and Ry is the 
hydraulic radius at any depth of flow Yo. By letting 


ACR2/3 =K 
2/3: 2 
AgCoRy?/3 =K, 


where K is called the conveyance of the channel, then S 


and S$ 
become : c 


f 
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Pees 
f K2 
Si 2 G: P 
Ot 2 
Ky 
d so 
Be So 
a hs (215 
ext, let o 
Be Or 8 


Pubstituting Fqs. 21 and 22 into’Fq. 20 gives 


2 
(1-4) (2&)¢ 
ee dy (23) 
K 
1- —%)? 
K 


‘quation 23 is a useful form of Fq. 13 and is used as the basis 
or discussing the various methods of integrating the equation 
f varied flow. 


a. Graphical Integration. 


This may be called an exact method, since no simplifying 
assumptions are involved. Equation 23 may be written in the form 
Sod ee d 
= (6! + — 
ee )? - 1 : 
50 (24) 
Integrating Eq. 24 between two sections x, and x, at which the 
depths are y, and y, gives 


il YQ 1-G 
Pert? axl f dy] 
ey So 1 Hy (Rye ] 
Kg 
or 
1 Yo 1] =6 
Mn a= Xo-Xy = a [ly = ¥,) + - pare bel (25) 
: 1 G)2-1 
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2 

To evaluate the integral in Eq. 25, the reciprocal ona - 
eee ; P 
is plotted as a function of (1-G)y in Fig. 5. 

For any known*values of y, and y, the shaded area in Fig. § 
obtained by a planimeter or analytically, gives the value of thi 
integral and the distance Ax may be calculated. By assumin' 
various successive values of y, the required curve is obtaine 
to any degree of accuracy desired. 


/ x 
(/-—G)y _— 


Fig. 5. Solution of Gradually Varied Flow Equation by Graphical Integration 


b. The Step Method 


This method of integration is widely used and is the simples; 
to perform. Equation 13 may be written 


6 ,Q? 
(ei herr eRe 
Ax = —————— Ay (26) 
So mS 


2 
; : 0) 
Provided Ay is small, an approximation of (1-28) and S, ca 
rs ‘ 


be obtained by evaluating these quantities as the mean of theii 
values at the two extremes of Ax. In other words, if ¥y, 1s the 
depth at x,, and y, is the depth at x, + Ax, then these quanti- 
ties can be evaluated by using a value of y = (y) + ¥9)/2. Singe 
Ay = yo - yy, the right side of Eq. 26 is completely known, ané 
the length of the reach Ax can be determined. By taking success- 
ive increments of Ay, the complete surface profile may be plot- 
ted. This method gives surprisingly accurate results. 
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¢.. Methods of Direct Integration 


The Method of Bakhmeteff.‘°’ Bakhmeteff developed a method for 
the integration of the varied-flow equation based on the use 
of a "varied-flow function" derived from two assumptions: 


(1) _ The conveyance of a channel K = ACR2/3, may be approxi - 
mated within a limited range of depths by the relationship 


K? = A*c2R4/3 = a constant x yN (27a) 
so that 
Logit yl SUG ener (27b) 
a Ag*Cy7Ry*/3 Y9 
where N is called the "hydraulic exponent." 
(2) Ge as ett and changes very slowly in any practical 


case so that for a reasonably short range of integration, it may 
be approximated by an average constant value. 

Using these assumptions and setting y/y,) = r and dy = Yy ar, 
Fq. 24 can be expanded and integrated from a section where x = x, 
yY = y, to a section where x = x5, y = yy resulting in 


Yo ye ar 
Ax = ES Sat ee tour yi = tree) a ag ere 


(28) 


; 

d 

The integral (fiers designated as B(N,r), the Bakhmeteft 
=f 


form as 


Ae = xy ~ xy =22f(rg-ry) - (1-6) (BIN,r) -BIN,ry)}} (29) 


The values of the varied-flow function covering a wide range of 
values of N and r have been computed and presented in tabular 
form'?’ to facilitate practical profile computations. In order 
to determine the surface profile in any given case using Eq. 29, 

the following step procedure is used: 
(1) determine the class of curve by evaluating y, and yg 
(2) determine the value of the "hydraulic exponent" N by 
means of a logarithmic plot of Eq. 27a. 
(3) plot the value of (1 - G) against y so that the average 
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value of (1 - G) between the assumed depths of flow 
may be evaluated 

(4) obtain the proper value of B(N, r) from the tables for: 
the terminal values of r 

(5) substitute the required values in Eq. 29 and compute 
the values -of.Az. 


The Method of Mononobe. ‘"? Mononobe developed a method for 
integrating the varied-flow equation based on the assumption that 
both the area A and the wetted perimeter p, may be expressed as 
monomial functions of the depth y. That is: 


Pp, = a constant x y! (30) 


A? = a constant x y@ (31) 


In Eq. 30 and 31 the exponents 1 and & will generally var 
slightly with the depth y; in these cases, however, their average 
values in the range of depths under consideration may be used., 
Mononobe has shown that the above empirical relationships are 
applicable to many types of cross-sections in common use. Using 
these relationships and Manning’s friction-slope formula (Eq 
llb) it can be shown that 


S y .4/31-5/3M 
— = (—-) 


0 Yo 


By substituting this relationship and other relationships derived 
from the assumptions in Eq. 23, Mononobe developed the followin 
equation of gradually varied flow: 
£0 Cp-m-1l 

V,2M " 
“ dr - k's v2 5 dr] (32) 
1 Cos $e teers A IE | 


y 
Ax = Xor-X) = ae mr 


in which r = y/y,, and Cy) = 5/3M - 4/31. Mononobe further sim-- 
plified Eq. 30 by using 


r ab 
r & 
Cy-M-1 
r r 
ers 


These are referred to as Mononobe’s "backwater functions." Upon 
substitution, Eq. 30 reduces to 


Bul. 404. FLOW IN MILD, UNIFORM CHANNELS 29 


Ax = x9-x) = —*{$,(r,) - O)(r,) -2 [4,(r,)-8,(r,)]} (35) 


The value of 6, depends only on the values of Comendien, 
whereas the value of &, 1s a function of Cy, M andr. The values 
of the backwater functions have been computed for various values 
‘of Cy, M and r and are available in graphical form.'’’ These 
charts are for use only in the case of backwater curves since 
the values of ) and ®, are given for corresponding values of r 
greater than unity. 

For convenience in integration in the cases of the dropdown 


curve, the variable is taken to be — instead of r. Accordingly, 
Eq. 32 is modified to ss 


1 -2 ,] M-1 
(S (=) 
ie ie 4 2y l/ro r 1 
eee Of 2 gt Caper) ie 8) 
So l/r, a eeu FH aa 1 ony ie 


fe 


and Mononobe’s dropdown functions are defined by 


=9 
Ara fg merle 
Pie ee e001 oer (37) 
r 
l 1 M-1 
ZS GF) 1 
W feed te ee eee 38 
27 ee etc Ti io 
ae) eee! 
4 that Eq. 36 may be written as 
Oy 2 
eee ee Oy eee ye 1) (89) 


As is the case with the backwater functions, the dropdown 
function y, depends only on the value of C, and 7; whereas fy 1s 
a function of Cy, M and>> The values of the dropdown functions 


corresponding to various values of C), M and i are given by 
Mononobe in graphical form. ‘+®? Z 
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In applying Mononobe’s method to profile computations it will. 
be found that logarithmic plots of the wetted perimeter against. 
depth do not approximate straight lines very well; therefore, an. 
error in the results may be introduced by assuming the wetted! 
perimeter to be a monomial function of the depth. The greatest, 
drawback in this method, however, lies in the inherent error due: 
to interpolation and the difficulty in using the charts to find. 
the values of the dropdown and the backwater functions. It will 
be noted from Eqs. 35 and 39 that the values of these functions 
are multiplied by the term y,/S) in computing the distance along; 
the channel Ax. The value of y)/S) is generally large, so that a. 
slight error due to interpolation may result in an appreciable 
error in the computed distance along the channel. For this 
reason, it has been suggested that the values of the varied-flow’ 
functions be given to the fourth decimal place in order to keep 
the error due to interpolation to a minimum. The graphical 
charts prepared by Mononobe are drawn on such a scale that the 
values of the backwater and dropdown functions can be read only 
to the second decimal place. It would appear that unless the 
charts are drawn to a much larger scale, the presentation of the 
values of these functions in tabular form is to be preferred. 


INTEGRATION OF THE GRADUALLY VARIED-FLOW EQUATION: 
PROPOSED METHOD 


- Assumptions 


In the foregoing discussion it is pointed out that assumptions 
€ necessary to simplify the solution of the gradually varied- 
w equation by direct integration. The equation is inte- 
ated, in each case, between the range of depths in which the 
implifying assumptions are valid. Thus, the number of steps 
quired in the integration procedure for the whole range of 
pths under consideration will depend on the limiting conditions 
f each set of assumptions. For instance in the "step method," 
e assumption that the mean friction slope is equal to the mean 
f the friction slopes at the ends of a reach is valid only for a 
all increment of depth. Accordingly, a comparatively large 
ber of steps is required in the computation procedure. In the 
ethod of Bakhmeteff the empirical relationship between the con- 
eyance and depth holds for a wide range of depths; therefore, 
or those cases involving insignificant effect due to velocity 
ead changes, the computations may generally be made in one step. 
f the effect of velocity head change is not negligible on the 
ther hand, the limitation imposed by assuming G to be constant 
akes it necessary to confine the computations to small incre- 
ents of depth. It is evident that the assumptions made must 
onform to the actual conditions over a wide range of limits if 
he method were to yield reasonably accurate results without 
ecourse to a large number of intermediate steps. 

In the light of the studies made by Bakhmeteff'?’ and Mono- 
obe,‘’’ it appears that both the area and the conveyance can be 
xpressed as monomial functions of the depth over a considerable 
ange for the channel sections in common use. These assumptions 
re made in the proposed method of integrating the gradually 
aried-flow equation, namely 


K2 = a constant x y% (27a) 


A2 = a constant x y/ (31) 


quations 27a and 3] are only approximate; some cross-sections 
atisfy these empirical relationships more closely than others. 


31 
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However, even in channels not well adapted to the monomi 
approximations over their entire range, the surface curv 
usually cover only a limited range so that the actual variati 
of the hydraulic elements conforms closely to these empirical re 
lationships. In most cases it will be possible to choose value 
of the constants and the exponents that will make Eqs. 27a and 3 
express the values of K and A with sufficient accuracy over 
wide range of depths. If greater accuracy in results is desired 
the whole range of depths under consideration may be subdivid 
to determine the values of the exponents separately for ea 
interval. However, this procedure will hardly be required ij 
making practical profile computations. 

It will be noted that the values of the exponents N and M ar 
characteristics only of the given cross-section. As such, th 
can be determined once and for all for a given channel. Hence, 
saving of time as well as labor of computations may be effect 
in those problems in which a large number of surface curv 
are to be computed for the same channel under a wide variatio: 
of discharge. 


13. Derivation of Proposed Equation 


The differential equation of gradually varied flow may 
solved with the aid of the simplifying assumptions by starti 
from Eq. 13% 1seé., 


bQ? 2 aA 2 
(1-285) oe apie yore 
dx = g dy & ay i 2g A 
So-S; Sy-S; SoS, 
ROP gies 
dy 2g Op 


. de—) as alqn) 

avd ee = + — 

rae 1-t (40) 
oe cr 
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erefore, Fq. 40 becomes 
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nd integration yields 


A 2 
dts?) V2 dG) 
So x=—+ {/——+ ie + constant 
Pome oo (Nal -28¥0, 2 (Fp yN 


33 


(41) 


(42a) 


(42b) 
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bere? & 
Designating as before y/y = r and B(N,r) = Iz=49r Eq. 42? 
becomes ia 


A.2 
d (9) 
S V2 A 
OR eat B(N,r) + 0 ——— + constant (43) 
Yo BY 9 1- (20) 
J 


The integral on the right side of Eq. 43 may be reduced to th 
form of the Bakhmeteff varied-flow function with the aid of Eq. 31 
from which it can be deduced that 


(Yoyw 2 ¢Xoyayayn, Aoynym 
hi ¥ A2 


Therefore A, 2 ee 
tle: ees 
Fishery ee (44) 
TON | y= ON ie 
y A2 


Now, the integral on the right side of Eq. 44 is in exactl 
the same form as Bakhmeteff’s varied-flow function as define 
previously. In this case the exponent is N/M and the variable 1i 


Ansa? therefore 


2 
SO 9 
A2 NA 
ee ia 
A,2 M’ A2 (45) 
1-(9-)N/m 
2 


and Eq. 43 becomes 

V 2 2 

x oo Ot BI Nord + ro ope osm) (46) | 
So 28Y 6 M A2 


in which the constant of integration has been intentionall 
omitted because Eq. 46 is used only to determine the diffe tan 
in the values of x between any two given sections. 

For the downstream section at which r = ro, and A = Ao, 


Fe vey ye iG A, (47a) 
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nd for the upstream section at which r = r,;, and A = A 


Lee 
rail (47b) 


2 
Vo 


28Y 0 


0 
My Ess [ry -B(N,r 


+ 


N 
) i 
’ BL, 


sing the A notation and subtracting gives 
2 


V N 
Ax = Xo~Xy Seer = ABU rt 420 Api, 


2 
Ao 


is )] (48) 


In applying Eq. 48 to practical profile computations, it is 
ecessary to know the values of both of these varied-flow func- 
Fions. Bakhmeteff has computed the first function for various 
yalues of r in the range of N lying between 2.8 and 5.4.‘ For 
he common types of practical channel sections the value of the 
pxponent N/M generally lies between 1.0 and 2.0 which is less 
han the range covered by Bakhmeteff’s tables. Hence a new set 
f tables has to be computed for evaluating the second varied- 
low function in this range of N/M, and for appropriate values of 
the variable Age Ae. This has been done with the results pre- 
ented in Tables 2 and 3 in which the values of the function are 
iven for sufficiently small intervals of both the variable and 
he exponent to permit interpolation. The values of the varied- 
low function for the range of Ag iA between 0 and 0.990 are 
iven in Table 2, and for the range between 1.01 and 20 in 
able 3. In Table 3 all the values shown have been increased by 

to prevent the values from changing signs. This procedure is 
ermissible since the relative values of the functions, which are 
sed in the computations, remain unaffected. 

The proposed method, as well as that of Mononobe, differs from 
akhmeteff’s method essentially in the treatment of the velocity 
ead changes. In Bakhmeteff’s method the factor G, reflecting 
hese changes, is assumed to be constant. The limits of integra- 
ion should, therefore, be confined to small increments of depth 
inasmuch as the factor G is a variable with respect to the depth. 
n the proposed method G is, in effect, regarded as a function of 
he depth so that the integration limits may be generally taken 
over the entire range of depths considered. For these cases, the 
proposed method will afford a more direct computation procedure 
whereby the results can be obtained without recourse to succes- 
sive steps. 

The integral embodying the effect of changes of velocity head 
in the proposed method (Eq. 45) is reduced to the form of the 
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varied-flow function defined by Bakhmetefl. This form is preferred 
in order to avoid the necessary adoption of a new term, as 1 

done in the method of Mononobe. In fact, both Mononobe’s back- 
water and dropdown functions are reducible to the same form as 
Bakhmeteff’s varied-flow function. It seems that little, wf 
anything, is gained by defining them in their present form, 
especially when the resulting integrals ®, and W,. are expressed 
in terms of three parameters. The integral in Eq. 45, on the 
other hand, is dependent upon only two parameters, N/M andg 


Table 2 
NA? Ad? 
Values of the Varied Flow Function B(—,-$-): —<| 
MA A 
N/M 
Ay?/A? 
1,00 1.20 1.33 1.50 1.67 2.00 
0.000 0.000 0.000 0.000 0.000 0.000 0.000 
0.05 0.051 0.050 0.050 0.050 0.050 0.050 
0.10 0.105 0.103 0.102 0.101 0.100 0.100 
OF25 0.163 0.159 0.156 0.153 0.151 O-15F 
0.20 0.223 0.216 0.212 0.207 0.203 0.203 
0.25 0.288 0.275 0.269 0.263 0.259 0.255 
0.30 0.357 0.337 0.329 0.322 0.318 0.310 
0.35 0.431 0.407 0.395 0.383 0.375 0.365 
0.40 0.511 0.480 0.464 0.448 0.437 0.424 
0.45 0.598 0.556 0.536 O.SLF 0.504 0.485 
0.50 0.693 0.642 0.617 0.592 0.574 0.549 
0.54 (ORCC 0.714 0.684 0.655 0.634 0.604 
0.58 0.868 0.804 0.759 0.725 0.700 0.663 
0.62 0.968 0.884 0.842 0.799 0.767 0.725 
0.66 1.079 0.979 0.930 0.881 0.845 0.793 
0.70 1, 204 1.086 1.029 0.973 0.931 0.867 
O43 1,309 bea Irae 1.109 1.049 1.005 0.929 
0.76 1.427 1, 281 1. 208 1,134 1.078 0.996 
0.79 1.561 1.387 1.305 1, 226 1, 168 1.071 
0.82 Lats 13525 1.430 1,333 1.260 a ep Sy 
: 
0.85 1,897 1.670 1.566 1, 467 1, 393 1, 256 
0.88 2.120 1,857 iis 1.619 13532 1.376 
0.90 2.303 2.019 1,880 1.741 1. 636 1.472 
0.91 2.408 2.103 1.959 1.817 1,710 1,528 
0.92 2.526 2.201 2.046 1.894 1.780 1.589 
0.93 2.659 2.314 2.148 1.984 1.860 1,658 
0.94 2.813 2.469 2.285 2.088 1,934 1.738 
0.95 2.996 2.598 2.405 Zc2l2 2. 066 1.832 
0.96 3.219 2.824 2.604 2.361 2.169 1.946 
0.97 Salt Oye 3.039 2.801 2.554 2.364 2.092 
0.98 3.912 3, 323 3.066 2.832 2.661 2.298 
0.985 4, 200 3.610 3.318 3.021 2.793 2.443 
0.99 4.605 4.044 3.699 3,294 2.966 2.647 
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9° /A?; consequently, Tables 2 and 3 representing the values of 
this integral are more simple to use than the graphical charts 
prepared by Mononobe'’?’ for the , and wy. functions. 


4i4. Computation Procedure 


In applying Eq. 48 to profile computations, it should be borne 
in mind that the downstream conditions constitute the upper limit 
of integration. The computations should be made with strict 
observation to sign. If the location of the depths is assumed 
correctly for the upper and lower limits, the distance between 
any two sections Ax will always appear as positive. * 

To compute the surface profile of flow in a given channel 
When the rate of discharge is known, the following procedure 
may be used: 

a. Determine, from the given conditions, the points of 
control from which the computations may start, and the type of 
surface curve that will develop. 

b. Determine the values of the exponents N and M in Eqs. 27a 
and 31 for the entire range of depths covered by the surface 
curve. This can be done algebraically by substituting into these 
equations a pair of simultaneous values of the depths at the ends 
of the range under consideration and the corresponding values of 
the conveyance or the area. It can also be done graphically by 
plotting on logarithmic scale the conveyance and area as ordi- 
mates and the corresponding values of the depth as abscissas. 
A straight line is then drawn to follow the points as closely as 
possible. The slope of this line is one-half the value of the 
exponent. The algebraic method corresponds to fitting a straight 
line through the initial and final points on the logarithmic 
plot; it has the disadvantage of not revealing the degree of 
approximation for the intermediate points. For this reason, the 
graphical method is more likely to give the average values of the 
exponents by following all the points as closely as possible. 

c. Compute the distance between any two given sections 
according to Eq. 48. The necessary parameters in finding the 

A 
values of the two varied-flow functions are (N, r), and Gras). 


The value of N will generally be between 2.8 and 4.0 and for this 
range the tables by Bakhmeteff‘?’ should be consulted. For the 
range of the exponent N/M between 1.0 and 2.0, Table 2 or 3 may 
be used, depending on the value of the variable. In the event 


*While either X) or both Xo and %, may have negative signs, Ax = X%y - %4 will 


always be positive. 


Table 3” 


N A? 
Values of the Varied Flow Function BO): Feat | 


Ps Dye! 

Ao a 1.00 1. 20 1.33 1.50 1.67 2.00 
1.01 7.605 6.944 6.621 6.295 6.046 5.652 
1,02 6.912 6.367 6.102 5.835 5.632 5.308 
1.03 6.507 6.033 5.801 5.566 5.387 5.107 
1.04 6.219 5.793 5.585 5.376 Dead 4.966 
1,05 5.996 5.609 5.420 5.229 5.084 4.857 
1.06 5.813 5. 456 5. 283 5.109 4,977 4.768 
1.07 5.659 5.328 5. 169 5.008 4.886 4.694 
1,08 5.526 jo 2a) 5.070 4.921 4.808 4.629 
1.09 5.408 ig pee 4.983 4.843 4.737 4,573 
1,10 5.303 5.045 4.913 4.775 4.669 4.522 
LG) 4.897 4.699 4,607 4,512 4,441 4.331 
1.20 4.609 4.462 4.394 4,329 4.282 4,199 
1525 4. 386 4, 286 4.236 4.187 4.151 4.099 
1.30 4.204 4,143 4.109 4.073 4.046 4.019 
1.40 3.916 3.901 3.897 3.896 3.896 3.896 
1,50 3.693 3.724 3.741 3.761 3.778 3.805 
1.60 3.01] 3.579 3.615 3.654 3. 687 3.733 
Ht) 3.357 3.457 3.509 3.564 3. 608 3.675 
1.80 3.223 3.358 3.423 3.488 3.538 3.626 
1.90 3.105 3.263 3.341 3.422 3.485 3.585 
2.00 3. 000 3.182 Sede 3.364 3.436 3.549 
2.10 2.905 3.108 3.209 3.312 3.392 3.518 
2.20 2.818 3.041 3.152 3.265 3.353 3.491 
2.30 2.738 2.983 3.102 3.223 3.316 3. 466 
2.40 2.663 2.928 3.056 3.185 3.283 3.444 
2.50 2.594 2.875 3.012 3.149 3.254 3.424 
2.60 2.530 2.826 2.971 3.117 3.229 3.405 
2.70 2.469 2.782 2.934 3.087 3.204 3.389 
2.80 2.412 2.738 2.898 3.059 3.181 3.374 
2.90 2.358 2.700 2.866 3.032 3.158 3.360 
3.00 2.307 2.662 2.835 3.008 3.139 3.347 
3.10 2,258 2.626 2.805 2.985 3.121 3.335 
3.20 Zeek 2.593 2.778 2.963 3.103 3.323 
3.30 2.167 2.560 2.751 2.942 3.086 3.313 
3.40 2.124 2.531 SBP! 2.923 3.071 3.303 
3.50 2.084 2.503 2.704 2.904 3.054 3.294 
3.60 2.044 2.471 2.679 2.887 3.043 3.285 
3.70 2.007 2.455 2.658 2.870 3.029 3.277 
3.80 1.970 2.419 2.637 2.854 3.017 3.269 
3.90 1.935 2.397 2.619 2.839 3.003 3. 262 
4.00 1.901 2.372 2.599 2.824 2.993 3.255 
4.10 1.869 2.348 2.580 2.810 2.982 3.249 
4.20 1.837 2.327 2.563 2.797 2.972 3.243 
4.30 1. 806 2.307 2.547 2.784 "2.960 3.237 
4.40 1.776 2.285 2.529 2.771 2.951 3.231 


“All values in this table have been increased by ‘3 to prevent the numbers from being 
negative. 


1.00 


1,747 
1,719 
1.692 
1.665 
1.639 


1.614 
1.565 
1.518 
1.474 
1.431 


1.391 
1.351 
1.314 
67400 
1.242 


1.208 
fhe Ate) 
1.144 
1.113 
1.083 


1.054 
1.026 
0.998 
0.972 
0.946 


0.921 
0.860 
0.803 
0.749 
0.697 
0.649 


0.602 
0.558 
0.515 
0.474 
0.435 
0.397 


0. 361 
0.326 
0.292 
0.259 
0.227 
0.197 


0. 167 
0.110 
0.056 


Table 3, (concluded) 


N/M 


eH 


33 


.514 
. 500 
. 484 
- 469 
- 455 


Nnwwnd bd 


- 442 
-416 
392 
370 
347 


NNN hd ds 


. 326 
. 308 
287 
. 269 
253 


236 
221 
205 
SOM 
176 


NnmNwwhd NONNWN ND 
° Spiers . . 


164 
149 
137 
~125 
113 


101 
074 
050 
025 
005 
983 


FPNmnwND bdo NNNNND 
. Oe is . CRSA 


965 
945 
927 
912 
. 896 
881 


a 


. 866 


oll eal oll eel anil eee 
amine’, 6 
co co OC 
eSNoPo 
ANrIoCp 


. 804 


= 
4 
Ko) 
_ 


1.773 
1.753 


ry 
uo 
So 


. 760 
. 748 
2130 
. 126 
LG 


. 106 
. 687 
- 669 
652 
- 636 


621 
607 
.594 
581 
. 569 


NOnNnwNw bd NNNN DY mONNWD 


5 
. 546 
535 
2020 
~515 


NNNNnNN 


506 
497 
. 488 
479 
471 


NnwWNN Nd 


- 463 
~ 445 
428 
412 
. 398 
2.384 


NNNMNW Nd 
* 


2.371 
Zoo 
2.348 
2.338 
2.328 
2.318 


2.309 
2.301 
2.292 
2.285 
PA PH 
2.270 


2.263 
2.251 
2.239 


NNNNwN 
. . e. . e 
a | 
~ 
_ 


— 
lon 
~ 


NNONMNwN bd NnwnwNWwh 
a Mispiie ia Cake he See ie 
co 2) 

_ iw) 

w nn 


NOnwnwhd wy 
2, eles eo) vie) 
co 
— 

[ae} 


NOnNnwNwNH hd Nw hd 
suite ntreiy ache ye ae tet 8m 8) 
4 
i) 

i} 


NN NYNNNNNDY NNNNNND 
; sorte nike reams ils! Sa snesen commie 
a 
= 
wn 


2.00 


3.226 
3.221 
3.216 
Seal 
3.207 


3.203 
3.195 
3.187 
3.180 
3.174 


3.168 
3.163 
3.158 
3.153 
3.148 


3.144 
3.140 
3.136 
3. 132 
3.129 


3.126 
3.123 
3.120 
Sai 
3.114 


Sek 
3. 106 
3.100 
3.096 
3.091 
3.087 


3.083 
3.080 
3.077 
3.074 
3.072 
3.069 


3.067 
3.065 
3.063 
3.061 
3.059 
3.057 


3.056 
3.053 
3.050 


40 ILLINOIS ENGINEERING EXPERIMENT STATION 


that either the exponent or the variable lies between the value 

shown in these tables, the value of the varied-flow function ma 

be found by straight-line interpolation. 
This computation procedure is illustrated by the following 


examples. 


a. Example 1. 

Let it be desired to find the distance along a rectangula 
10-ft-wide channel between the depths y, = 7.0 ft and y, = 4.5 f 
when the rate of discharge is 136.0 cu ft per sec. The slope o 
the channel S, is 0.0004, and the value of n, the roughnes 
factor in Manning’s equation, may be taken as 0.015. 

~ The uniform depth, yg, computed from Manning’s equation 11 
4.0 ft for the given rate of discharge. The critical depth y, 11 
computed to be 1.79 ft. The profile of flow, therefore, belong 
to the case of backwater, or the M,-type of surface curve. 


Table 4 
Hydraulic Properties of the 10-ft Rectangular Channel 


Wetted Hydraulic 
Depth Area Perimeter Radius Conveyance 
y A Py R R2/3 K 
2.0 20 14.0 1.43 1.28 2,540 
3.0 30 16.0 1.87 ioe 4,530 
4.0 40 18.0 2.22 Lgl 6, 800 
5.0 50 20.0 2.50 1.85 9,190 
6.0 60 22.0 ate 1.96 11, 630 
7.0 70 24.0 2.92 2.05 14, 200 
8.0 80 26.0 3.08 16,900 


The hydraulic properties of the 10-ft rectangular channel ar 
shown in Table 4, in which the conveyance K is computed accordin 
to Eq. lla. The value of the hydraulic exponent N is found fro 
the logarithmic plot of conveyance-depth shown in Fig. 6. The 
slope of the line is 1.40; therefore the hydraulic exponent N has: 
a value of 2,80. The value of M in Eq. 31 is 2.0, since the area 
is directly proportional to the depth for channels of rectangular 
cross-section, . 

The given data are 


Geer) <2 ROR i 
So 0.0004 10; 000 
ep? > Sch hy 8 , Reo 


28Yo 64.4 x 4 
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erefore 0 
Bree eT oe— ot ya" -0.625 
om the tables prepared by Bakhmeteff, 


BE2;6,.16125)..= 0.705 
B(2.8, 1.75) Qe 212 


nce 
AB(N,r) = -0.492 


e parameters necessary for using Table 2 are 


20000 


a7 


6 000 


Conveyarice “FT 
A 
N 
S 
S 


£000 


1000 7 


a BEY i A AES GLY) 


3 
Depth a7 9 


Fig. 6. Logarithmic Plot of Conveyance K Against Depth y 
for the 10-ft Rectangular Channel 


4I 
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M 25.0 
2 
Heal Renter) 
AS? 45° 
A, * 40? 
— = = 0.33 
2 
me 70 
from which 
BCLEA, 0979) ale 273 
B(1.4,0.53)-= 0.365 
Therefore 9 
A 


Neo 
ABGT*75-) = -0.908 


Substituting these values into Eq. 48 gives 


Ax = x. - x, = 10,000 [0.625 + 0.493 - 0.045 (0.908)] 
= 10,770 ft 


b. Example 2. 


Determine the distance along the rectangular channel in tk 
previous example between the depth y, = 2.0 ft and y, = 3.6 ft 
The profile of flow in this case obviously belongs to the dropdown 
or Mj-type of surface. The given data are 


Therefore 


Ar @irg = ry, 994040 
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om Bakhmeteff’s tables, 
B(2.8,0.90) = 1.253 
B(2.8,0.50) = 0.521 


i] 


nce 


AB(N,r) = -0.732 


ye parameters necessary for using Table 3 are 


A.2 402 

— = — = 1,24 
Ay? 362 

A.2 402 
aoe 
An? 20 

ee a 1 AO 


om which 
DeICA 1524). =A, 246 
6(1.4,4.00) = 2.692 


herefore 


ubstituting these values into Eq. 48 gives 


Ax = x. - x, = 10,000 [-0.40 + 0.732 - 0.045 (1.554)] 


2620 ft 


. Example 3, 


Let it be desired to determine the extent of the dropdown 
irve in a 10-ft circular conduit (n = 0.015) laid on a slope of 
001 when the rate of flow is 305 cu ft per sec. 
imit of the curve may be taken as 0.95 of the uniform depth. 
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The upper 


The uniform and critical depths are first computed for the 
iven rate of flow and are found to be 6.0 ft and 4.11 ft respec- 


ively. The values of the exponents N and M may be determined 


scording to Eqs. 27a and 31 as follows: 
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Taking logarithms of both sides of the equations yields 


K A 
loge) logy) 
oer Capea moras Ee 
l logG— 
os) og) 


The exponents may be found by substituting into these equations 
pair of simultaneous values of the depth at the ends of the rang 
under consideration and the corresponding values of the convey 


ance or area. 
In the example above when y = 0.95 x 6.0 = 5.7 ft 


K = 8900; A = 46.20 
and when yy, = 4.11 ft 


Ky = 5060; Ay = 30.42 


Therefore 


108900 


and 


The required distance along the channel may be determined: 
according to the scheme of computations in the following table., 
Applying the results in the table to Eq. 48 yields, 


Az = zy - x, = (S0, [(0,685 - 0.95) - (0.733 - 1.369) 


2 
+ ae (2.986 - 4,65)] 


= 1230 ft 
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Computations for Dropdown Curve in 10-ft Circular Channel 


Upstream Section, yioe oat LE Downstream Section, Yo =°4. 11ete 
y 
f= 22 ='0,95 rae) 22) & 0,685 
Yo Yo 
B(N,r) = 1.369 B(N, ry) = 0.733 
A,2 (49.2)? A,2 2 
Qe = 1.14 pte 422) = 2.63 
Ay? (46.2)? Boome (c0.42)" 
BR A” = 4.650 Be Ay" 2.9 
ran : Bra? = . 986 
. Example 4. 


Supposing the 10-ft circular channel in the previous example 
S preceded by a chute in which the flow of 305 cu ft per sec 
akes place at a uniform depth of 2.0 ft, what will be the length 
f the surface curve along the channel? 

The depth of flow at the junction between the circular channel 
nd the chute will be taken as 2.0 ft. Since the flow must 
cross the critical depth of 4.11 ft to attain a uniform depth of 
5.0 ft in the circular conduit, the surface curve will end in a 
1ydraulic jump. The problem, then, is to determine the length of 
the M,-type surface curve that lies between the limiting depths 
of 2.0 ft and the lower conjugate depth of 2.7 ft at the foot of 
the jump. The computations are shown in the following table, 
assuming the same values for N and M as in Example 3. 

Substituting these results into Eq. 48 gives 


Computations for the Length of the Mg~-Type Surface Curve 
in the 10-ft Circular Channel 


Upstream Section, y, = 2.0 ft Ay (49.2)2 
ary = <arer = 19.60 
A (11.18) 
rune | = 0.333 ; 
Yo 2 
ab ae 1.846 
B(N,r}) a 0.334 MA 2 _ ; 
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2 2 
Downstream Section, y> = 2.7. ft Ag, (0:2)se piers 
A, (18 Ve 
lo = Y2 = 0.45 
Yo N Ag? 
BO—>) = 2.205 
B(N, ro) = 0.456 M A,? 
: BoA, 2 : 2 
Sixt wea yet ae 0.001 [(0.45 - 0.333) (0.456 - 0.334) 
+ (6.2)? (2,205 - 1.846)] 
64.4% 6 
= 186 ft 


The application of the proposed method to profile computatio 
has been illustrated by solving examples involving the thr 
types of surface curves of flow in uniform channels on mil 
slopes. In all the foregoing computations the velocity distri 
bution is assumed to be unity, and no account is taken of t 
eddy losses accompanying retarded flow as in the case of t 
backwater and M3-type of surface profile. Equation 48 may bt 
modified to include the effect of nonuniform velocity distributio 
as well as the eddy losses, which are usually estimated a: 
some percentage of the change of velocity head between adjacens 
sections. The modified form of the equation, following thi 
procedure in deriving Eq. 48, becomes 


Yo Vy, tas uae 4 : 

Ax =— [Or = ABIN yr) 4 (CG e-"k) ABE ae (49) 
0 287 

in which the correction factor for eddy losses k and the momentur 
factor C, are assumed to be constant from section to sectio 
It should be noted that these factors tend to offset each othe 
due to their opposite sign; the net effect, generally, is t 
prolong the length of the M, curves and shorten the M-type o 
surface curves. For instance, using values of C. = 1.05 an 
k = 0.10 (assumed) in the solution of Example 1 would yield 
result of 10,870 ft instead of 10,770 ft; and in Example 4, th 
length of the M, curve would be 174 ft instead of 186 ft. Thes 
results show only a small deviation from those computed by mean 
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of Eq. 48. In the event, however, that the corrections for non- 
uniform velocity distribution and eddy losses are large enough 
to be significant, the computations may be made in accordance 


with Eq. 49. 


15. Special Charts 


In the foregoing examples the distance along the channel is 
computed between two sections of known depths for a given rate of 
flow. In other types of problems it may be required to determine 
|the depth of flow at the end of a given reach of channel. When 
Ax or y is to be computed over a wide variation of the discharge, 
much of the subsequent labor of computation can be avoided by 
preparing an auxiliary curve for the solution of the gradually 
varied-flow equation. This can be done by plotting x“0 against 
r according to Eq. 46* Y0 


Sy V52 A,2 
*2— = r= BIN, rn) - pee eye lanp ee 
Yo 8Y 9 M A 

Vo7 
a rf N .-M) 
r B(N,r) Deyo oye” 


In the above equation the values of the exponents N and M are 
‘determined by the given channel section; therefore, it is possi- 


| ; : 
ble to compute the values of x SO corresponding to the various 


2 YO 
0 3 A 
values of ogy -and r for the channel under consideration. 
| 0 
=) 
Such a plot is shown in Figs. 8a and 10a, in which x 2 is 
Y0 


plotted against the variable r and a parameter C, for various 


values of the exponents N and M. It is shown in the following 
discussion that the particular values of the exponents used for 
these figures correspond, respectively, to the special cases of 
channels with rectangular and circular cross-sections. 


*From Eq. 31 


An? 
eee eS imi ro 
A Yo 
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a, Rectangular Channels. 

For channels of rectangular cross-sections the area is direct- 
ly proportional to the depth; therefore, the exponent M has a. 
value of 2.0 according to Eq. 31. The hydraulic exponent N may’ 
be determined in accordance with Fqs. lla and 27a. 


K = 1.486 4R?2/3 = a constant x y"/2 
n 


Thus it can be seen that the value of N can be determined either: 
by a logarithmic plot of the conveyance against depth as ex- 


plained previously, or of the product AR 2/3 against depth as; 
shown in Fig. 9. In the latter case, the roughness factor n need 
not be taken into consideration in computing N. 


Table 5 
Values of the Exponents M and N for Rectangular 
Channels of Various Widths 


Hydraulic R Dpepth 
Width, b Radius,.R a aE cae M N 

ft ft From To 
b<y Red ee So eS Cee 
5.0 if 1.0 8.0 2.0 2.60 

b 

10.0 R= Aree 2.0 12.0 2.0 2.70 
20.0 4.0 14.0 2.0 2.80 
30.0 4.0 16.04 | ..2:0 » 2) aye0 
b>>y R®y 2.0 3.33 


For rectangular channels of broad cross-section the hydraulic 
radius R is roughly proportional to the depth. Therefore, the 


product AR2/3 varies approximately as the 5/3 power of the depth, 
in which case N is equal to 3.33. The minimum value of N is 
attained in rectangular channels whose width is small compared 
to the depth, for which the hydraulic radius remains relatively 
constant. Consequently, N is equal to 2.0 in the case of rect- 
angular channels with deep and narrow cross-section. For the 
intermediate cases N lies between these limiting values. 

Values of N determined from a logarithmic plot of the product 
AR2/3 against depth y for the various widths of channel are shown 
in Table 5. The average value of N is found in each case for the 
range of depths stated. 
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From the results in Table 5 it can be seen that the exponents 

N and M for a broad rectangular channel are, respectively, 3.33 

and 2.0. Using these particular values of the exponents, the 
S 

values of x—™ calculated from Eq. 46 are shown in Table 6 for 
0 

various values of r >. 1 and for the parameter C,. The results 


are also presented graphically in Fig. 7, in which the value of 
So. 
Bet separated into two parts, [r - B(N,r)] and CBee ree) 
0 


for convenience in plotting. 


rae S 
Similarly the values of x —2 for N and M equal to 2.8 and 2.0, 

Y0 
corresponding to an average width channel, are computed in Table 
¢ and plotted in Fig. 8. The values for r >. 1 are in Fig. 8a 


where again the components [r - B(N,r)] and CBEE eM) are used 


for convenience. In Fig. 8b are the values for r <1. The 
application of Fig. 8a and 8b is recommended for rectangular 
channels between 10 and 30 ft wide since the hydraulic exponent 
in this range deviates only slightly from the average value of 
2.8 used in the computations. 


6. Circular Channels. 


In the case of circular channels the exponents M and N may be 
found by a logarithmic plot of the area and the product AR2/3, 
respectively, against the depth. This has been done in Fig. 9 
for 5, 10, and 20-ft circular conduits in a range of depths 
between 0.2 and 0.7 of the diameter. The plot shows the slope of 
the lines to be remarkably consistent with one another, the ex- 
ponent N having an average value of 3.60 and M, 2.56. Applying 
these values of the exponents to Eq. 46, it is possible to con- 
struct the auxiliary curve for solving the equation of gradually 
varied flow in circular channels. Table 8 shows the computations 


of ak for various values of the parameter C, and variable r. 
Yo 
The results for r >.1 are plotted in Fig. 10a in which the value 


of ieee again separated into the two parts. Figure 10b shows 
Y0 

the curves for r < 1. The curves in Fig. 10a and 10b may be 

used to determine the surface profiles of flow in circular channels 

between 5.0 and 20 ft in diameter. 
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In applying the charts to profile determinations it should be: 
noted that the curves possess several characteristic features, 
First, all the variables appear as ratios; the resulting plot is; 
therefore dimensionless. Second, the use of the curves 1s not: 
confined to any one channel roughness since the exponents N and M/ 
are determined only by the shape of the channel section. Third, 
the curves for r < 1.0 (Figs. 8b, 10b) possess a maximum point 
corresponding to the value of r, at the critical depth. Since: 
this point constitutes the lower limit of the dropdown curve, the? 
portion of the auxiliary curve comprising values of r greater? 
than r, is to be used in connection the the M,-type surface) 


curve; the remaining part of the curve, less than r,, is to be: 
used for the M,-type surface curve. In the event that a given 


Ses : 
value of x—° is satisfied by more than one value of r, the: 


Y0 
required result can be obtained by considering the physical. 
conditions of flow. 


i6. Practical Application of the Charts 


The special charts will be found to effect a saving of time: 
and labor of computations in those problems involving a wide: 
variation of discharge, as shown in the following examples. 


a. Example 5 


A natural channel, whose cross-section may be taken as a broad 
rectangle, is laid on an average slope of 0.0004. Uniform flow 
is attained at a depth of 4.0 ft when the discharge is found to) 
be 12.0 cu ft per sec per unit width of channel. Find the dis-: 
tance between two sections, 6 and 9 ft deep respectively, wheni 
the discharge is 15 cu ft per sec per unit width, assuming: (a) 
the effect of nonuniform velocity distribution and eddy losses to) 
be negligible; (b) the effect of velocity head changes to be 
negligible; and (c) c, = 1.10, k = 0.20. 

The uniform depth for q = 15.0 cu ft per sec may be found 
as follows: 


Q q K So’ y 38 


t 


ae err a | 

Q% 9 Kg So Yo 
Therefore 

3/5 3/5 
q bo 

Yo = yt) = 465) = 4,87 fe 
Hence 

V2 ren.) 1 q? 

Coe ee eG 
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Fld ae 


Graphical Solution of the Equation of Gradually Varied Flow 
in Broad Rectangular Channels 
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ee pose 


« oa =r-8(3.6,r)+CB(14, r7*) 


(N=360, /1=256, r=/) 
em 


Fig. 8. Graphical Solution of the Equation of 
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The distance along. the channel may be found with the aid c 
Fig. 8. The parameters required in using the chart are r and C; 
In Example 5(a), 

V.2 


0 
C, = Dgy, ~0.0365 


pan ge Pl ee 


Yo 4.57 


Applying these values of the parameters to Fig. 8 yields 
r - B(N,r) = 1.030 


-M 
c,BY,r ) = 0.025 
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Fig. 9. Logarithmic Plot of the Area A and of the Product AR2/3 
Against Depth y for Circular Channels 
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AC wp) 


B(2.8,°)+CB(.4,r4) 


(N=2.8, M 
rire 


—= r— 


So 
Yo 


x 


(4'EDG-4 


10. Graphical Solution of the Equation of 


Fig. 
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Therefore i 
0 
17, 1:030 + 0-025 =Fi035 
Y2 
Similarly, at the downstream section where fy = 4900 1.9% 
0 
S 
x, —2 = 1.878 + 0.008 = 1.886 
Yo 
Whence 
S S S 
Ax —& = xy — - x —= 1.886 - 1.055 = 0.831 
Y0 Yo Y0 
or 


Axo 20110. 831) =) 9500 14 
So 


For the case of negligible velocity head changes, the integral 
; . -M : : 
embodying the kinetic energy changes; i.e., BCE r ), 1S omitte 


from the final equation (Eq. 48). In Example 5b, therefore, 
Cc, =0 and 


Anca. TAr =. ABUSE 34a rl 
So 


11,430 (1.878 - 1.030) = 9690 ft 


For Example 5c in which the effects due to non-uniform velo 
city distribution and eddy losses are not negligible, the modi- 
fied equation of gradually varied flow should be used. Referrin 
to Eq. 49 it can be seen that 


V.2 
C, = (C, - k) 2— = (1.10 - 0.20) (0.0365) = 0.033 
28Yo 


At the upstream section where r, = 1,31, 
r=" BUN, r)).=" 1,030 


BBtreria) = 0; 028 


At the downstream section where ry Mrlioe 


r - BIN,r) = 1.878 


C,BUE, r=) = 0.0075 
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Therefore 
a 
Ax = — (1.878 + 0.0075 - 1.030 - 0.023) = 9,520 ft 
0 
b. Example 6 


Suppose the rate of flow in the broad rectangular channel in 
the previous example is to be measured by two gaging stations 
5000 ft apart. If the stage at the upstream gage is 5.0 ft and 
the fall is 1.5 ft as indicated by the downstream gage, what is 
the rate of flow in the channel? Assume C= 1.10 and k = 0.20. 

The problem consists of selecting a discharge such that the 
fall between the two sections will be equal to 1.5 ft' under the 
given conditions. This can be done by trial, or by graphical 
interpolation from a plot of the fall against assumed rates of 
low as follows: 

The maximum fall between the two sections is attained when 
the surface is level, for which 


Ce) 
=9¥,.* 2000'S, 


SS 
I<) 
i 


Ngan yy 2 2.0 ine 
The minimum fall occurs where the flow is uniform, 1.e., 


OF Shy 55. 0 rit 
. 9/3 

fae C5 

: 4 
Since the fall of 1.5 ft lies between the maximum and minimum 
conditions, the corresponding flow must assume an intermediate 
value between zero and 17.4 cu ft per sec. To plot the fall- 
discharge curve, the computations may proceed by taking a series 


= Lied tefs 


Computation of the Fall-Discharge 


ft assumed Yo Yo Yo ft ft 
cfs/ft 
Q. = Se eee level surface - - - ---- - 7.00 2.00 


sie O,029 iG NO) OE LPO A PHE SK — Lhatehil 
: BS 02030 1.43" 15236" 07571 T2807] 1.90 6:65 “1.65 
12.00 4.0 0.032--1.25" 0.951 “0.500 1.451 1.60 6.405 1.40 
4.5 01033. 1. kl) 05623 0.444_ 1.067 16327 —5.94" 0,98 

AAO)! 2 a ee uniform flow - - ------- 5.00 0 
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2.0 
Fige Il fFall-Discharge Curve 
in Example 6 
D/.0 
q \ 
25 /0 £0 


va 


of intermediate flows and finding the corresponding fall, the 
upstream depth y, being held constant at 5.0 ft. This is done ir 
the preceding table. 

From these results it is apparent that the flow corresponding 
to a fall of 1.5 ft lies between 9,60 and 12.0 cu ft per sec. By 
graphical interpolation from the fall-discharge curve in Fig. ll 
it is found to be 11.5 cu ft per sec. 

If it 1s desired to establish the stage fall-discharge re- 
lationship for the channel similar curves may be drawn, eack 
corresponding to a fixed depth at the upstream section. Fro 
such a plot it is possible to determine the discharge in the 


channel when the upstream gage reading, together with the fall, 
1S given. 


c. Example 7. 


Water is fed from an impounding reservoir to a storage reser- 
voir by gravity flow through the rectangular channel described ir 
Example 1. The length of the channel Ax is 5000 ft, the slope 
Sy is 0.0004. Assuming the depth at the outlet of the impounding 
reservoir to stand relatively constant at 5.0 ft, construct the 
delivery curve of the channel for various depths at the inlet tc 
the storage reservoir. 

The delivery curve (see Fig. 12a) has several special points 
At a there is no flow, since the water surface is level. Thi 
corresponding depth at the lower section is 


Yo = 5.0 + (0.0004) 5000 = 7.0 ft 


At point 6b, corresponding to y, = 5.0 ft, the flow is uniform ana 
the rate of discharge is 


: 
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Q = KS* = 9190 x 0.02 = 183.8 cfs (from Table 4) 


. point c, corresponding to Yo = ¥,, the delivery is a maximum; 
9 additional flow can be obtained by further lowering y, below 
1e critical depth y,- Therefore, for the portion of the curve 
ying between points a and b, corresponding to rates of flow 
slow 183.8 cu ft per sec, the surface profile belongs to the 
,-type for which Fig. 8a is applicable. At higher rates of 
yw the surface profile assumes the M, curve, for which Fig. 8b 
ay be used. 

The computation procedure is illustrated by first determining 
particular point on the delivery curve, lying between points a 
id 6. For this portion of the curve the rate of discharge is 
sss than 183.8 cu ft per sec. 

Assume Q = 136 cu ft per sec, for which 


Yo =o4.() ft 
V 2 
C, = ~*~ = 0.045 
28Y 9 
5 
ra sat oe ly Oe 
Yo 4 


S d ’ 
1e value of x —°, read from Fig. 8a, is 0.793. Therefore 
Y9 


U0 
x, = 0.793 ea 7930 ft 
sequently 
oan ne + x, = 5000 + 7930 = 12,930 ft 


» that, 
201 12,930 (0.0001) = 1.293 
0 
om Fig. 8a, ry, is found to be 1.55. Hence, the depth at the 


wer section is 


yo = 1.55 (4.0) = 6.20 ft 


a 


For a point on the delivery curve between points 6 and c, the 
ite of flow is greater than 183.8 cu ft per sec. Assume Q = 204 


ft per sec, for which 
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0 SO 100 [$0 @00 &50 
Discharge Q 
(6) 


Fig. 12a. Delivery of a Canal; Fig. I12b. Delivery of 
the 10-ft Rectangular Channel in Example 7 


Yo = 9.4 ft 
V 2 
c, = — = 0.041 
28Y 
pip egt ls weet §2095 
Yo 
S 
The corresponding value of x—® from Fig. 8b is -0.26. Therefore, 
Y0 
= -0.26 (2:4 _) = -3510 ft 
“1 F:0008 
Consequently, 


%o = Ax + x, = 5000 - 3510 = 1490 ft 
and, 


S 
0 0.0004 

= 1490 (¥: = 0.11 
*2 5 ere 


From Fig. 8b it is seen that for C, = 0.041, the maximum value 


S ‘ 
that x —2 may reach is 0.093, corresponding to r, = 0.433, them 
Yo RS 
being no value or r which would satisfy x—2 = 0.11, The physi- 
Y0 


cal meaning of this result is that the M, curve corresponding to 
Yo = 9.4 ft and Q = 204 cu ft per sec is shorter than the lengtk 


of the channel. In other words, the delivery Q = 204 cu ft 


per sec is in excess of Qe. 


To find the value of the maximum delivery, which obviously 
lies between 204 cu ft per sec and the uniform flow of 183.8 cu ff 
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er sec, it 1s necessary to determine the values of vy, for 
everal intermediate rates of flow, making in each case, Yo = VG: 
e discharge, which results in ¥, = 5.0 ft is the maximum de- 


livery. The method of computation is shown in the following 
table, frow which the maximum discharge is found to be 202 cu ft 
per sec. ; 


Computation of Maximum Delivery 


Ss Sj S, 
Q Yo Cy ve ro oe 0 peel ry Yi 
assumed Yo Yo Yo 
cfs ft ft ft 
188 OF O42 e222 ee ORAS5 0.094 0.393 -0.299 0.936 4.77 


5.1 
194 a2 0.042 2.27 0:5436 0.094 0.385 -0.291 0.935 4.86 
paso We04l)  *2.34 © 0.436 05,092 =70.574 © 0-282 07932585;00 
The computations of y. for different rates of flow in the range 


between the maximum delivery and Q = 70.0 cu ft per sec are 
tabulated below. Assuming y, to remain constant at 5.0 ft, the 


delivery of the channel for various depths y, at the lower 
section is plotted in Fig. 12b. 


Computation of the Delivery 


(y; = 5.0 ft) 
S § S 

Q sie Yo Cy ry eee pa) pee lo ¥y 
cfs ft ft Y0 Yo Yo ft 
0 -+--+-+-+-+-+---- level surface - - ----+-+-+----- 7.00 
70.0 tet t 2e5 0.049 2. 00 215045 0.800 Dey PRU  Waisiil 
oA .5 108 2. 23..,55 0.046 ah) ae Oswell 1685 aCe Ontos 
136.0 Le) Caw 0.045 2 Sa Oe 93 0.500 12293, led SamOscU 
160.0 2.00 4.5 0.044 A (a2 Pal 0.445 0,872 61:29 26-80 
183.8 2.20 5.0 - - - - uniform flow ------------- 5.00 
188.0 ED RSE | 0.042 0.980 -0.644 Oe393°% .-0.251. 0.9250 4am 
194.0 LAC Bee 0.042 0.961 -0.440 O¢385 .-0.055. 0.695 3.6] 
202.0 9.34 5.36 -- - - maximum flow ------------- Pay: 


The delivery curve in this example is obtained under the 
assumption that the depth at the upper section y, remains fixed. 
In the general case, however, both the depths y, and y, will 


vary. The problem, then, is to determine the rate of flow in the 
channel for any possible combination of depths at the given 
sections. This can be done by computing as before a series of 
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delivery curves, each corresponding to a fixed value of the depth 
at one of the given sections. The final results may be repre- 
sented graphically on a three-dimensional plot, showing the 
relationship between the fixed depth, fall and discharge. 


d. Example 8. 
The flow in a 10-ft circular conduit (S,) = 0.001, n = 0.015) 


is 305 cu ft per sec. The depth at a certain section is found to: 
be 5.70 ft. Compute the depth at a section 1230 ft downstream 
from the first section. 

This problem is obviously the reverse of Example 3. The: 
uniform and critical depths are respectively y, = 6.0 ft and 
y. = 4.11 ft. The other data are 


c, = = 0.10 
ey mR 
= 2:1. 0,95 
Pe. 6 a eines 
: ; S 
Using these results in Fig. 10b, x —2 is found to be 0.05. 
Therefore Y0 


Yo 


—° (0.05) = 300 ft 
So 


ro ee 


x» = 300 + 1230 = 1530 ft 


and 


S 
xo—2 = 1530 (0-001) = 9,255 
Yo 6 


From Fig. 10b it is found that for C, = 0.10, x So ~~ 0. 255.am 
satisfied by two values of r; viz, 0.72 and 0.63. The correspon- 
ding depths are, respectively, 4.32 ft and 3.78 ft. Evidentl: 
the required depth is 4,32 ft since the dropdown curve cannot 
cross the critical depth of 4,11 ft to reach a depth of 3.78 ft 
without passing through a hydraulic drop, a condition to whic 
the gradually varied-flow equation no longer applies. 

It should be noted that these results are obtained with the 
aid of Fig. 10b, which is based on N = 3.60 and M = 2.56. If the 
actual values of the exponents for the 10-ft circular channel are. 
used (N : 3.50, M = 2.58), the required depth should come out the 
same as in Example 3; viz, 4.11 ft. The deviation of results ini 
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using the prepared chart is roughly 5.0 percent, which falls 


within the degree of accuracy warranted by most practical profile 
computations. 


17. Summary 


The application of the proposed method to profile determina- 
tions has been clearly illustrated by solving several typical 
examples. Special charts based on the proposed method were drawn 
to facilitate the solution of the equation of gradually varied 
flow in rectangular and circular channels. These charts are best 
adapted to the determination of the delivery in a canal, thereby 
affecting a great saving of time, as well as the labor of compu- 
tation, over the graphical or the step methods. The computa- 
tions, from which the special charts were prepared, are also 
given. Using these computations the charts should be reproduced 
to a more usable scale. In the event that the actual exponents 
for the given channel vary slightly from those used in preparing 
Figs. 7, 8, and 10, the charts may be applied to yield at least a 
first approximation, as shown in Example 8. For channels of 
other cross-sections similar charts should be prepared according 
to the particular values of the exponents. 


Vl. LABORATORY INVESTIGATIONS 


18. Purpose and Scope 


The purpose of the laboratory investigation was twofold: to 
determine the deviation of the actual surface profiles from 
the theoretical; and to obtain data on actual surface profiles; 
as the basis for comparison between the proposed method andl 
the various methods of computation that are in common use. 

In the laboratory the actual surface profiles of flow in ex-: 
perimental channels of rectangular, trapezoidal and circular: 
cross-sections were recorded for different rates of flow. The: 
channels were all set on a sustaining slope of 0.001 due to the: 
limited capacities of the pumps used. Only the backwater (M,-: 


type surface curve) and the dropdown (M,-type surface curve) were: 


observed, inasmuch as they are the more common types of surface: 
curves encountered in practice. The length of the backwater: 
curve was limited by the depth for uniform flow and the height of! 
the side walls of the experimental channels; therefore, a low 
rate of flow was used in the backwater experiments to provide ai} 
wider range of depths covered by the M, curves. Similarly, ai 
high rate of flow was used in the dropdown experiments since, in! 
this case, the length of the M, curve is confined to the range: 
between the critical depth and that for uniform flow. For this; 
reason, each surface curve was observed only for the rate of flow/ 
most suitable under the limitations of the experimental equip-: 
ment. The different rates of discharge used in the experiments. 
are given in Tables 10-12 and ranged between 0.176 cu ft per sec 
and 0.425 cu ft per sec. 


19. Apparatus 


In the arrangement shown in Fig. 13 two 2-in. centrifugal. 
pumps were used to lift water from the sump into the constant-: 
head tank (A) from which the water flowed by gravity into the: 
‘weir box (B) through a 4-in. connecting pipe. The rate of flow 
was controlled by means of a 4-in. gate valve, the excess water: 
being returned to the sump from the constant-head tank through a 
3-in. drain. The water flowed over a 90-deg V-notch weir in the 
weir box, and thence fell freely into the entrance tank (C). 
hanging baffle was provided in the entrance tank to cut down the 
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inlet turbulence. The depth of flow in the experimental channel 
(D) was controlled by adjusting the opening of the tail-wate 
gate in the outlet box (£), from which the water returned t 
the sump. 

Two motor-driven centrifugal pumps were used. The combinec 
capacity of. the pumps was approximately 200 gal per min under ¢ 
20-ft head. The two 2-in. discharge pipes from these pumps met 
in a 4-in. riser pipe connected to the bottom of the constant- 
head tank. 

The constant-head tank consisted of a 50-gallon drum placec 
concentrically inside a larger cylindrical tank, 30 in. ir 
diameter. The rim of the drum, approximately 6.0 ft around, 
served to maintain practically a constant head for the rates oj 
flow used in the experiments. The excess water overflowed the 
constant-head tank into the large tank and returned to the sump 
through a 3-in. drain. To regulate the rate of flow, a 4-1in. 
gate valve was installed in the pipe connecting the constant-heaa 
tank to the weir box. 

The weir box was constructed as a rectangular tank 2 x 2.5 x» 
ft. A gravel-filled baffle was placed at the inlet end of the 
weir box to reduce the entrance turbulence of the incoming flow, 
The tank was separated into two compartments by a 90-deg V-notck 
weir which was used for measuring the rate of discharge. The 
location of the weir plate, made of 1/8-in. black iron, and the 
dimensions of the V-notch weir are ‘shown in Fig. 13. The rate o) 
discharge corresponding to a given head on the weir was read fro 
a calibration curve. The head on the V-notch was measured by : 
point gage set above the stilling well which was connected to thi 
side of the weir box at a point 14 in. upstream from the weir: 
A 14-in. rectangular opening extending the full depth of the wei. 
was cut in the downstream side to permit a free fall of wate: 
into the entrance tank. 

The entrance tank consisted of a 30-in. drum, in which 
diametrical wooden partition extending 2.5 ft into the tank wa 
fastened perpendicular to the direction of flow. The partitio) 
served as a hanging baffle to reduce the inlet turbulence as th 
water entered the experimental channel. 

The depth of flow in the channel was controlled by adjustin 
the tail-water gate in the outlet box which was mounted at th 
downstream end of the channel. The tail-water gate consisted o 
three vertical vanes made of wood with short pieces of shee 
metal projecting from the edges as guide vanes. The wooden vanea 
were pivoted at the bottom of the outlet box. By turning t 
vanes about their vertical axes the area of the opening w 
adjusted to produce the desired depths of flow in the channel 
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view of the experimental channel and the outlet box with the 
aneS 1n position is shown in Fige 4. 

Point gages were used to measure the depths of flow along the 
hannel and the head on the 90-deg V-notch weir. These gages 
Pig. 15) were made of 3/8-in. round steel rods, 15 in. long, and 
hreaded for a distance of 12 in. The unthreaded portion was 
apered to a point; the top of the threaded portion was shaped 
ad fitted with a small circular handle for turning the gage. 
1e side of the threaded stem was flattened to record the number 


fthe threads from the handle. The base of the gage, a 3/4 1] 


Fig. I4. End-View of Channel Showing the Tail-Water Gate 
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Fig. 15. Point Gage 


x 6-in. block of steel, was tapped and fastened rigidly to 
wooden frame across the top of the channel. The face of the bas 
was marked with ten equal divisions along the diameter of th 
rod. By counting the number of threads on the rod above the bas 
and the fraction of a thread on the base-scale, it was possibl 
to take readings to the nearest one-tenth of a thread, that is,, 
to 0.006 in. The depth of water was found by taking the dif- 
ference of the gage readings at the water surface and som 
other reference point. 

The experimental channel itself was located above the sump, 
the elevation and plan of which are also shown in Fig. 13. I 
rested on a framework of two 2 x 6-in. beams, spaced 5 in. apar 
with 2 x 4-in. spacers. Side braces were used to hold th 
channel sections rigidly in place. The beams were supporte 
every 8 ft by 4 x 4-in. posts, which were braced laterally and 
longitudinally with 1 x 6-in. planks. A series of holes wer 
drilled in the posts through which a 1/2-in. round iron pipe was 
inserted cross-wise between each set of supports. A 4 x 4-in., 
crosspiece was placed on top of the iron pipe. The channel ang 
the supporting beams rested on sets of screw jacks, which in turn 
rested on the cross-pieces. The method of support is shown int 
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ig. 13. By moving the iron pipe vertically along the series of 
oles in the posts, the channel could be set on any desired 
lope. The final adjustments were made with the jack screws 
hich had a maximum lift of 2 in. A side view of one of the 
xperimental channels and the supports is given in Fig. 16. 

Three types of channel sections - rectangular, trapezoidal 
nd circular - were used in the tests. The rectangular channel 
as made of 14-gage black iron, bent into 8 x 8-in. sides; the 
Tapezoidal section was made from the same channel by bending the 
ertical sides 70 deg with the horizontal. These channels were 
l|@abricated in lengths of 4 ft and 8 ft. The joints of the 
eparate sections were made water-tight by caulking them with 
lasticine clay. The total length of the channels was 40 ft; 
Owever, to avoid the error caused by turbulence at the inlet 
nd outlet of the channel, only the’ 24-ft center section was used 
n observing the surface curves in the tests. 

The circular channel was made of 24-gage black iron in 10-ft 
ections. An additional framework was used to support the 
ircular channel. It consisted of rectangular boards, spaced 
% {t apart, each of which had a 10-in. diameter semicircle cut 
ut of the top. Rectangular slots, 1 x % in., were cut at the 
uarter points of the semicircle and longitudinal braces were in- 
erted through these slots along the entire length of the channel. 
he framework rested on the 2-in. beams, with spacers fitted in 


Fig. 16. Side View of Experimental Channel and Supports 
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between to prevent any horizontal movement, The sheet iron we 
rolled by hand in the laboratory into the framework to form th 
circular channel. Because of the inadequacy of the method usec 
kinks appeared in several places along the channel. The fina 
dimensions of the channel were determined by blocking the ends o 
a horizontal section, into which known volumes of water wer 
added and the corresponding depths accurately determined by mear 
of the point gages. The average area of the cross-section at 
given depth was found by dividing the corresponding volume 

water by the length of the section. The results, as well as t 
actual values of the wetted perimeter, are given in Table 

The total length of the circular channel was 40 ft, of which onl 
the 22.5-ft center section was used as the effective length 11 
observing the surface profiles. 

The dimensions and the hydraulic properties of the experime 
tal channels are shown in Table 9. In determining the surfa 
roughness of the channels, the uniform depths correspondin 
to four different rates of discharge were observed in each cha 
nel, the results being tabulated below. The value of n, t 
roughness factor, 1s computed by substituting the given data int 
Manning’s equation. In making subsequent profile computation 
for the various channels the average value of n, given in thi 
last column of Table 9 is used. 


20. Experimental Procedure 


Each channel was set on a sustaining slope of 0.001. This w 
done by letting water stand in the entire channel and adjusti 
the jack screws until the depth of water, measured with the poi 
gages, increased uniformly in the downstream direction. Sin 
the point gages were spaced every 8 ft, the increment of dept 
between the gages was adjusted to 0.008 ft for the desired slop 

The depth of water in the channel was measured by lowering t 
gages until the point barely touched the surface. The differenc: 
between the gage reading at the water surface and the invert o 
the channel was converted to the depth by a conversion factort 
which was determined by comparison with a standard gage. Thi 
head on the weir was found in a similar manner, except that ii 
this case the zero reading was taken at the crest of the V-notc: 
weir. The rate of discharge used in the tests was found from 
calibration curve for the recorded head on the weir, which wa: 
checked before and after each run. 

The laboratory procedure used in determining the sur fac: 
curves followed essentially the same technique introduced b: 
Mononobe.'*°! That is, the entire surface curve was obtained by 
conducting a series of profile determinations, the initial dept! 
1n a new run being adjusted approximately to the final dept: 
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Table 9 
Hydraulic Properties of Experimental Channels 


——————————— a eae a 
————————————— ODO ewww 


Wetted Hydraulic 
= Depth, y Area, A Perimeter, Py Radius, R Conveyance, K G 
eG sq ft £t ft 


se a a ce 
Rectangular Channel 
Slope, So = 0.001; Roughness, n = 0.0092 
ee ee er eee se 


0.24 0.160 1.147 0.140 6.952 0.246 
0.28 0.187 227 0.152 8.597 0.234 
0.32 0.213 1,307 0.163 10. 273 0.228 
0.36 0. 240 1.387 Onis 12.041 0.219 
0.40 0.267 1. 467 0.182 13.861 0.211 
0.44 0.293 1,547 0.189 15.608 0. 200 
0.48 0.320 1.627 0.197 17.506 0.194 


4 a 
Trapezoidal Channel 


Slope, Sp = 0.001; Roughness, n = 0.0106 
a ee ae ee es ee 


0.24 0.181 3 lise Mate 0.154 7.287 0.235 
0.28 0.215 1.263 De ileal 9.298 0.235 
0.32 0.251 1.348 0.186 11.460 0.233 
0.36 0. 287 1.433 0.200 13.770 0.231 
0.40 0.325 1.518 0.214 16.279 0.229 
0.44 0.364 1.603 0.227 18.957 0.229 
0.48 0.404 1.688 0.239 21.794 0. 228 
0.52 0. 445 ds 0.251 24.834 0.228 
0.56 @. 488 1.858 0.262 27.966 0.223 
Circular Channel 
Slope, So = 0.001; Roughness, n = 0.0102 
0.099 0. 030 0.508 0.060 0.669 0.228 
0. 134 0.048 0.596 0.081 ol 0.254 
0.166 0.065 0.680 0.096 1.983 0.255 
0.202 0.088 0.773 0.114 3.008 0.260 
0.235 0.110 0.849 0.130 4.103 0.268 
0.268 0.132 0.919 0.144 Daze 0. 268 
0.299 0.155 0.987 0.157 6.573 0.273 
0.328 0.177 1.052 0.168 Ege 0.273 
0.357 0.199 1.109 0.179 9,214 0.273 
0. 387 Or22i 1.193 0.185 10.464 0.265 
0.414 0. 243 1.240 0.196 11.952 0.265 


in the run before. By connecting the series of profiles thus 
ybtained, it was possible to extend the total length of the sur- 
face curve beyond the limit confined to the effective length of 
che experimental channels. The applicability of this experimen- 
-al procedure has been well substantiated by tests conducted by 
fitchell and Barron. ‘13! 

In the profile tests the flow was first adjusted to uniform 
conditions by manipulating the tail-water gate. The opening of 
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the gate was then changed by making a small turn of the vanes t 
establish nonuniform flow in the channel. The surface profil 
was recorded by four point-gages, spaced 8 ft apart for the rect 
angular and trapezoidal channels and 7.5 ft apart for the cire 

lar channels. The readings were taken only after the flow has 
become steady. In the succeeding runs the upstream depth wai 
adjusted to duplicate the downstream depth in the run be fores 
This procedure was repeated until the flow reached a deptt 
slightly smaller than the height of the channel in the backwateg 
experiments, and larger than the critical depth in the dropdo 

experiments. In the circular channel the flow was too turbulen 
to permit accurate readings of the gages at the smaller depthss 


Table 10 


Comparison of Test Data in Rectangular Channel with 
Results Computed by Various Methods 


Distance Test Graphical Step Bakhmeteff Proposed 
Et. in. in. in. in. in. 


Backwater 


Sp = 9.001 Q= 0.176 cu ft/sec Yo = 2.44 in. 


0 3.02 3.02 ee 3.02 Sie Oy 
24 3.20 3.19 e198 aoe EAL 
48 sTou 3.36 3.36 3.36 3.36 
72 3.58 3.54 3.54 3.56 3. OG 
96 3.78 Bat Sig 0b Bark 3.78 

120 3.99 S297 3.98 4.00 4.01 
144 4.24 4.20 4.21 4.25 4.23 
168 4.45 4.44 4.42 4.41 4.48 
192 4.66 4.68 4.68 4,72 4.72 
216 4.91 4.93 4.91 4.97 . 

240 5.16 5.20 5.19 Deed a 
264 5.42 5.45 5.45 5.49 aS 
288 5.70 Sd ale 5.76 Salt 

: 
Dropdown 


So = 0.001 Q = 0.346 cu ft/sec Jou 359 ins 


3.84 3.84 3.84 3.84 
3.80 3.80 3.82 3.80 
3340 3.76 3.78 3.75 
3.69 3.69 3.75 Seite 
364 3.60 3. 68 3. 66 
3.50 3.50 3.59 Bao 
3.35 3.34 3.49 3.47 
3.08 3.03 3.35 3.33 
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[o reduce observational error the initial depth in the succeed- 
ing run was made slightly larger (dropdown) or smaller (back- 
water) than the final depth in the preceding run. 


21. Results and Discussion 


The results of the profile determinations in the various 
channels are presented an’ Tables 10, 11 and 1257 in’ which the 
abserved depths are given for the corresponding aie pe eee along 
she channel. The test data are also presented in graphical form 
oy plotting the elevation above an arbitrary datum against the 
distance. The datum, indicated in Fig. 17, is arbitrarily chosen 
at a distance 5 in. below the origin 0, Peter ne to the initial 
section from which the observed profiles are recorded. Therefore 
the elevation corresponding to a given depth y at a distance x 
from the origin may be computed according to 


El. = 5+ y - Sox 


Table || 


Comparison of Test Data in Trapezoidal Channel with 
Results Computed by Various Methods 


Test Graphical Step Bakhmeteff Proposed 
in. in. in. in. in. 


Distance 
it. 


Backwater 


Sp = 0.001 Q = 0.205 cu ft/sec Cae 2.74 in. 


0 3.42 3.42 3.42 3.42 3.42 
24 3.61 3. 60 3.60 3.58 3.59 
48 3.82 3.81 3.80 3200 3.78 
72 | 4.03 4,02 4.02 3.98 Sinks] 
96 4, 28 4.25 4,25 4,21 4,22 
120 4,49 4,49 4, 48 4.44 4,46 

144 4,72 4.73 4.73 4.70 4.70 
168 4.97 4.99 4.98 4.93 4.94 


Dropdown 


So = 0.001 Q = 0.425 cu ft/sec yo = 4.34 in. 


WwWwWwwWHw. hp 
wOwwwnw pp» Ps 


CWWww hap 
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The elevations, computed from the observed depths, are plottec 
against the corresponding distance from the origin; the data are 
shown as small circles in the diagrams. Backwater measurements 
in the rectangular, trapezoidal and circular channels are pre- 
sented in Figs. 18a,b - 20a,b; corresponding dropdown measure- 
ments in the various channels are shown in Figs. 18c,d - 20c,dl 
It must be emphasized that these figures are drawn to a muck 
distorted scale, so that a slight experimental error may result 
in a noticeable deviation from the general trend of the data. 

The results, obtained from the profile determinations in the 
various experimental channels, are subject to the following 
sources of error: 


(1) Observation. In the dropdown experiments the surface 0) 


Table 12 


Comparison of Test Data in Circular Channel with 
Results Computed by Various Methods 


Distance Test Graphical Step Bakhmete ff Proposed 
it: in. in. in. in. in. 


Backwater 


Spy = 0.001 Q = 0.183 cu ft/sec Yq = 3.34 in. 


Dropdown 
= 0.001 Q = 0.318 cu ft/sec yo = 4.53 in. : 


aonu-s 


. 
. 
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Cer > PPP RS 
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3. 


amoorn 
ne Ofo 


Bul. 404. FLOW IN MILD, UNIFORM CHANNELS 79 
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Fig. 17. Elevation above Datum 


le water became too turbulent to permit precise measurements of 
we depth. The fluctuations of the water surface were more pro- 
yunced in the circular channel. In addition, when the water 
wed past the joints of the channel sections’a rippling motion 
| the surface was produced which was carried some distance 
yond the joints. 

(2) Point Gages. The point gages were made in the laboratory. 
scause of the slight bending of the rods it was impossible 
) move them vertically without some side deflection of the 
int. However, since the side motion was limited, it is not 
lieved that the depth would be greatly affected as a result of 
ie slightly inclined positions of the rods. 

(3) Structural Defects. The circular channel was rolled by 
ind, giving rise to kinks in several places. The crown of the 
ivert was displaced from the longitudinal axis where the kinks 
‘re so formed that the location of the zero reading of some of 
ie gages could not be accurately determined. 

A slight deformation of the channel sections at the ends made 

difficult to align the joints to a smooth surface. A small 
ount of leakage took place at some of the clay-filled joints, 
t the losses were not appreciable compared to the rates of flow 
ed in the tests. 


After the channel was set to the desired slope, the deflection 
the section between beam supports was checked. This was done 
finding the depth at midspan below the level water surface. 
e results indicated that the assumption of uniform slope 
ulted in but a slight error. ; 
Theoretically, the observed data should lie on a continuous 
rve. It can be seen from Figs. 18-20 that the trend of the 
ual data was generally along a smooth curve. The deviation of 
e of the points, somewhat exaggerated in these figures due to 
distorted scale used, is probably caused primarily by the 
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observational error. The extent of experimental error, however 
is obviously fairly small in view of the close agreement betwee 
the observed and theoretical results as presented in Tables l 
disand tls 


22. Comparison of Actual and Theoretical Profiles 


The test data are compared with results computed by the graph 
ical, step, and Bakhmeteff methods, as well as the method prc 
posed in this bulletin. The theoretical curves were obtained | 
computing the distance along the channel between known sections 
In each case the initial section was taken at the origin 0 i 
Fig. 17, and the computations proceeded by small increments c 
depth. The results of the various methods of computation a1 
plotted in Figs. 18-20 for easy comparison with the test date 
The backwater curves in the various channels are shown separate} 
for the several methods. Those computed by the proposed ar, 
graphical methods are shown in Figs. 18a, 19a and 20a; those } 
the step and Bakhmeteff methods, in Figs. 18b, 19b and 20) 
Similarly the dropdown curves in the various channels for tl 
proposed and graphical methods are given in Figs. 18c, 19c an 
20c; those for the step and Bakhmeteff methods, in Figs. 18d, 19 
and 20d. The elevations above the arbitrary datum are re: 
directly from the backwater and dropdown curves for various dil 
tances along the channel. The results are also presented 
Tables 10, 11 and 12, in which the depths shown are computed fr 
the corresponding elevations according to Fig. 17. 

The method of Mononobe is omitted from these tables becau 
his charts'’) are not suitable for computing short sectio) 
such as those used in the observed profiles. For this purpo: 
the charts should be redrawn to such a scale as to be compatib 
with the required degree of accuracy of the results. This co 
not be done, because the original values of the varied- 
functions for the preparation of the charts are not inclu 
in Mononobe’s paper. 

Tables 10, 11 and 12 show that the discrepancy between 
observed profiles and the theoretical curves in the vari 
channels is reasonably small. The actual deviation is somewh 
exaggerated in Figs. 18-20 as a result of the distorted sc 
used in plotting the curves. In addition to experimental er 
the discrepancy between observed and theoretical profiles may 
due to any of the following causes: 


(1) Incorrect assumptions in computation procedure. An erz} 
may result from using average values of the exponents, N and 
and the factor G, in the proposed and Bakhmeteff methods; « 
from using the evaluation of the mean friction slope as t 
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average of the slopes at the ends of a section in the ste 
method. The error, however, is likely to be small because in t 
~ frst two methods the values of the exponents and G are relativ 
ly constant in the range of depths covered by the observ 
profiles (see Table 9) while in the third method a large numb 
of steps is used to give better accuracy. 

(2) Erroneous assumptions in derivation of varied-flow equa 
tion. In solving the differential equation of gradually vari 
flow by graphical integration no assumption is made in the integ 
ration procedure; therefore any discrepancy between the result 
of the graphical method and the test data, other than that due t 
experimental error, may be taken as evidence of the fallacy c 
the assumptions made in deriving the varied-flow equation. Thes 
assumptions are: 


(a) hydrostatic distribution of pressure intensity 


(b) uniform velocity distribution; i.e., C, = 1.0 


(c) losses confined mainly to friction, computed 
as if the flow were uniform under identical 
conditions of depth, discharge, and boundary 
roughness, 


Of these three assumptions the first 1s inconsequential excep 
in the immediate vicinity of the critical section. It seem 
unlikely that this assumption would lead to any appreciable erre 
inasmuch as the dropdown measurements were not extended to tk 
critical depth. The error in the results is probably due pri 
marily to the second two assumptions. 

The effects of the second assumption may be determined by 
consideration of the gradually varied-flow Eq. 13: 


b Q2 b O2 
Piao Pet oe bi 
d aA d EA d (50) 
x = ————— dy - ————-_ dy 
So - Sy So - Sy, 


The effect of taking velocity distribution into consideration i 
to increase the value of the second term on the right. It wil 
be noted that the first and second terms are always opposite i 
sign, the latter being negative, for the M,- and M,-types c 


surface profiles. Therefore the effect of nonuniform velocilt 
distribution is to shorten the distance along the channel betwee 
two given depths. It follows then that in the case of the bac 
water profile the actual curve will lie above the theoretic: 
one, resulting in depths larger than those computed by the grap 
ical method; while for the dropdown curve the actual profile wi. 
fall below the theoretical resulting in actual depths which ai 
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smaller than those computed by assuming uniform distribution o 
velocity (i.e., C, = 1.0). 

Finally, the assumption that the losses in gradually varie 
flow are the same as for uniform flow may lead to erroneous re 
sults. This is due to the fact that the flow in the laborator 
tests is subject to an additional loss from turbulence, caused 
the water flowing past the joints of the channel sections. Thi 
effect of this additional loss is to increase the friction slop; 
S, in Eq. 50, with the results that the difference between th 


two slopes, Sy and S,, becomes smaller in the case of the back 


water and larger in the case of the dropdown curves. Consequent 
ly, the distance along the channel between two given section 
will be increased for the backwater and reduced for the dropdow 
profiles. In either case the actual profile will lie below th 
theoretical curve, yielding depths smaller than those computed 
neglecting turbulent losses. 

Thus it can be seen that in the backwater case the effects 0: 
nonuniform velocity distribution and turbulence tend to offse 
each other. The net result will depend upon which of thes 
effects is predominant. For the backwater determination in th 
circular channel (Fig. 20a), the test data fall consistentl 
below the theoretical curve computed by the graphical method 
This is probably due to the losses caused by turbulence, a fact 
which is omitted in the theoretical profile computations. Fo 
the rectangular and trapezoidal channels (Figs. 18 and 19a), som 
of the actual data are greater than the theoretical values. Thi: 
discrepancy may be due to the effect of the nonuniform distribu 
tion of the velocity, which is neglected in the graphical metho 

In the dropdown case the effects both of nonuniform velocit 
distribution and turbulence tend to shorten the actual profil 
from the theoretical curve. Therefore the observed depth at 
given section should be smaller than that computed by the graphii 
cal method. This can be seen from the dropdown measurements ij 
the trapezoidal and circular channels, shown in Tables 1] and ] 
(or Figs. 19c,d and 20c,d). The greater discrepancy in result! 
for the dropdown curves may be caused, partly at least, by th 
effects of nonuniform distribution of velocity and turbulen 
losses, both of which produce errors of similar sign. 

It should be emphasized that the maximum deviation of the tes 
data from the theoretical values computed by the graphical methc 
is reasonably small, (of the order of 0.12 in.). This clo 
agreement may be taken as evidence that the assumptions in t 
gradually varied-flow equation did not result in any substantia 
discrepancy between the theoretical and the observed profiles. 
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The comparison between the experimental data and the result 
of the graphical, step, proposed, and Bakhmeteff methods, give 
in Tables 10, 11 and 12, show the degree of accuracy of t 
computations in that order. Both the graphical and the st 
methods are likely to yield more precise results because of t 
large number of intermediate steps taken in computing the surfac 
profiles. In the proposed and Bakhmeteff methods the results a 
not affected by the number of intermediate steps, which are us 
only to compute the elevations along the channel. In fact, thes 
methods tend to become less accurate when applied to sho 
reaches of channel, since in these cases a corresponding err 
due to interpolation of the varied-flow functions yields a wid 
discrepancy of results. The deviation in the proposed method 1 
probably caused by interpolation; however, thé error is reaso 
ably small in view of the close agreement between the actu 
profiles and those computed by the proposed method. In general 
there is only a slight difference between the results of t 
proposed and Bakhmeteff methods. This is due to the fact th 
the observed profiles are confined to a limited range of depth 
in which the actual variation of the factor Gis small. Und 
these conditions it makes little difference whether G is assum 
to be constant as in Bakhmeteff’s method, or considered as 
variable with depth as in the proposed method. 


11. CONCLUSIONS 


The accuracy of the proposed method is determined by comparing 
he computed surface curves with the observed profiles in various 
xperimental channels of rectangular, trapezoidal, and circular 
ross-sections. The test data are also compared with the results 
omputed by the other existing methods. The method of Mononobe 
S omitted because his graphical charts cannot be read with the 
esired degree of accuracy. From the results of the experimental 
Avestigations, the following conclusions are drawn: 


(1) No substantial error is introduced by the assumptions 
made in the gradually varied-flow equation as indicated 
by the close agreement between the observed profiles 
and those computed by the graphical method. 

(2) The surface curves computed by the proposed method 
agree reasonably well with the observed profiles so 
that it may be expected that the proposed |/method can be 
applied to practical computations with reasonable 
accuracy. 

(3) The computations of surface profiles by the graphical, 
step, proposed, and Bakhmeteff methods yield only 
slightly different results, the order of accuracy 
being the same as the order of the above list. 


It must be emphasized that the steps used in making the compu- 
ations were so small that the graphical and step methods yielded 
2sults which were more accurate than those obtained by the pro- 
ysed method. In practice, however, the proposed method will be 
ecidedly more convenient to use inasmuch as the computations 
sed not proceed by successive increments. The results of the 
xperiments indicate that the precision of the proposed method 
alls well within the degree of refinement warranted in most 
practical profile computations. 
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Vill. SUMMARY 


The basic assumptions leading to the equation of gradualll 
varied flow are discussed at some length. The simplification ii 
based upon the one-dimensional method of analysis, in which th 
flow is considered as being essentially parallel. The differ: 
ential equation for gradually varied flow is derived from th 
simplified equation of momentum, in contrast to the more common 
procedure of deriving it from energy considerations. The differ 
ence between the two forms of the equation is explained. It ii 
concluded that the form derived from the momentum principle it 
more consistent with the empirical expressions for the eval 
uation of the friction slope since the latter expressions ar 
essentially momentum equations. 

The solution of the differential equation for gradually varie 
flow by graphical, approximate and direct integration is discusse 
with regard to the advantages and limitations of each method 
The choice of the method to be used is decided by the type o 
problem under investigation. Theoretically, the graphical metho 
should yield results of the greatest accuracy inasmuch as n 
assumption 1s made in the integration procedure. The work in 
volved in this method, however, is quite tedious so that it ha 
not been widely used. In many practical problems the metho 
of approximate integration or the step method may be used wit: 
sufficient accuracy. The step method is limited by the fact tha 
computations must be made by small increments of depth in orde 
to satisfy the assumptions made in evaluating the mean frictio. 
slope within each interval. The objection is not serious, how 
ever, in solving the type of problems in which only a few surfac 
curves have to be computed and elevations are needed all alon 
the channel. In another type of problem, on the other hand, suc 
as the determination of the delivery in a channel connecting tw 
reservoirs, a large number of surface curves are required for 
larger number of discharge rates. Under these circumstances, t 
use of the direct integration methods will prove to be of decide 
advantage in the time and labor of computations saved. 
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ire incomplete in that either the effect due to velocity changes 
s neglected or the friction slope is evaluated by means of the 
hezy equation with a constant coefficient, the use of which 
as been outmoded. Moreover, these methods are based on some 
dealized channel cross-section of simple geometric shapes. 
hese simplifications may lead to a substantial deviation between 
he computed and observed profiles, as shown by the experiments 
onducted by Mononobe. For closer agreement between the theore- 
ical and actual results, the actual shape of the channel must be 
onsidered, and the friction slope should be evaluated by a more 
uitable friction formula than that of Chezy. . 

By using an empirical relationship between the conveyance 
md the depth, Bakhmeteff was able to derive a more general 
olution in which the effect due to the shape of the channel may 
taken into consideration, and the friction slope may be evalu- 
ited by any suitable friction formula. His method is convenient 
hen applied to the problems in which the effect of velocity head 
hhanges is insignificant. Otherwise, the correction that must be 
ade to take the velocity head changes into account necessarily 
lestroys the convenience of his method. In the method of Monono- 
e the computation procedure is more direct; unfortunately it 
s marred by the difficulty of using his graphical charts and 
he inherent approximations that are likely to result due to 
nterpolation of the varied-flow functions. 

These difficulties are avoided in the proposed method which 
S, in effect, an extension of Bakhmeteff’s method to cover those 
ases involving significant effects due to velocity head changes. 
he integral, embodying these changes, is reduced to the form of 
he varied-flow function defined by Bakhmeteff. The values of 
he integral have been computed over a wide range, sufficient to 
eet the conditions encountered in most practical profile compu- 
ations. The convenience of the proposed method is illustrated 
y solving several typical examples. 
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